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ABSTRACT 


Volumetric pattern analysis of fuselage-mounted airborne antennas 
at high frequencies is the object of this research. The primary goal 
of this investigation is to develop a theoretical solution for pre- 
dicting radiation patterns of airborne antennas in an accurate and 
efficient manner. This is an analytical study of airborne antenna 
pattern problems in which the antenna is mounted on the fuselage 
near the top or bottom. Since it is a study of general-type air- 
craft, the aircraft is modeled in its most basic form. The fuselage 
is assumed to be an infinitely long perfectly conducting elliptic 
cylinder in its cross-section and a composite elliptic cylinder in 
its elevation profile. The wing, cockpit, stabilizers (horizontal 
and vertical) and landing gear are modeled by "n" sided bent or 
flat plates which can be arbitrarily attached to the fuselage. 

The volumetric solution developed in this study utilizes two 
elliptic cylinders, namely, the roll plane and elevation plane 
models to approximate the principal surface profile (longitudinal and 
transverse) at the antenna location. With the belt concept and the 
aid of appropriate coordinate system transformations the solution 
can be used to predict the volumetric patterns of airborne antennas 
in an accurate and efficient manner. Applications of this solution 
to various airborne antenna problems show good agreement with scale 
model measurements. Extensive data are presented for a microwave 
landing antenna system. 


n 


CONTENTS 

Page 

LIST OF FIGURES iv 

Chapter 

I INTRODUCTION 1 

II THEORETICAL BACKGROUND 6 

A. Diffraction by a Wedge 7 

B. Diffraction by a Curved Surface 13 

C. Diffraction by an Infinitely 

Long Elliptic Cylinder 18 

D. Near Field Scattering by a 

Finite Bent Plate 27 

III MATHEMATICAL MODELING OF AIRCRAFT 

FOR PATTERN COMPUTATIONS 39 

A. Analysis of Elevation Plane Model 39 

B. Analysis of Roll Plane Model 52 

C. Three-Dimensional Model Approximation 

of an Aircraft 63 

IV VOLUMETRIC PATTERNS OF AIRBORNE ANTENNAS 98 

V MICROWAVE LANDING SYSTEM AIRBORNE ANTENNAS. 122 

VI SUMMARY AND CONCLUSIONS 164 

Appendix 

I COORDINATE SYSTEMS TRANSFORMATION. 167 

II GEOMETRICAL CONFIGURATION FOR OFF-CENTER LINE 

MOUNTED ANTENNAS.... 171 

REFERENCES 180 

i i i 



Page 


LIST OF FIGURES 

Figure 

1 Geometry for three-dimensional wedge 

diffraction problem i.... 9 

2 Two-dimensional wedge diffraction geometry 14 

3 Diffraction by a smooth curved surface. 15 

4 Diffraction by a smooth closed cylindrical surface... 19 

5 Geometry of antennas mounted on an infinitely 

long elliptic cylinder 21 

6 Diagram showing the elliptic cylinder 

coordinate system. 26 

7 Bent plate geometry.. 28 

8 (a) Single diffraction of incident field. 

(b) Single reflection of incident field... 32 

(c) Double reflection of incident field 

(d) Single diffraction of reflected field 33 

9 Eo radiation pattern for a small dipole mounted 
above a square plate for 9=90° and 0j<(j)£360° at 

f=l 0.436 GHz 34 

10 (a) H-plane pattern 35 

(b) E-plane pattern.... 36 

11 (a) H-plane pattern 37 

(b) E-plane pattern 38 

12 Theoretical model of KC-135 aircraft.,., 40 


13 Illustration of surface waves propagating across 

the junction line of a composite elliptic cylinder.,.. 42 

14 (a) Elevation plane pattern for a x/4 monopole 

mounted forward of the wings on a KC-135 
aircraft. (No radome and vertical 
stabilizer included). 

(b) Elevation plane pattern for a x/4 monopOle 
mounted above the wings on a KC-135 aircraft. 

(No radome and vertical stabilizer included) 46 


iy 



Page 


Figure 

15 (a) Elevation plane pattern for a circumferential 

KA-band waveguide mounted forward of the wing 
on a KC-135 aircraft. (No radome and vertical 
stabilizer included). 

(b) Elevation plane pattern for a circumferential 
KA-band waveguide mounted above the wings on 
a KC-135 aircraft. (No radome and vertical 
stabilizer included) 


16 (a) Elevation plane pattern for an axial KA-band 

waveguide mounted forward of the wings on a 
KC-135 aircraft. (No radome and vertical 
stabilizer included). 

(b) Elevation plane pattern for an axial KA-band 
waveguide mounted above the wings on a KC-135 
aircraft. (No radome and vertical 
stabilizer included) 


47 


48 


17 Approximation of an open-end waveguide by an 
array of 5 infinitesimal broad wall elements with 
different weights and an array of 3 infinitesimal 

narrow wall elements with uniform weight 49 

/ 

18 Truncated composite ellipse as simulation of 

fuselage with radome included 51 

19 Total field distribution 51 

20 (a) Scale model (1/35) of Space Shuttle in Anechoic 

Chamber at NASA ;(Hamptpn, Va.) 53 

(b) Theoretical model of ,1/35th scale model of 

space shuttle 54 

21 Elevation pattern of monopole on a l/35th scale 

model of space shuttle (with radome included) 55 

22 Elevation pattern of circumferential slot on a 

l/35th scale model of space shuttle (with radome 
included) 56 


23 Elevation pattern of axial slot on a l/38th 
scale model of space shuttle (with radome 

included) 57 

24 Elevation pattern of axial slot on a l/35th scale 

model of space shuttle (with radome being 2") 58 


V 



Figure Page 

25 (a) Reflection problem in X'-y plane. 

(b) Diffraction problem in y-z plane 59 

26 Fuselage and wing geometry for theoretical aircraft 
model looking from the front. The antenna is 

always mounted on top of the models... 61 

27 (a) Roll plane pattern (Ea) for a 1/25 scale model 

of a KC-135 with a A/4 monopole on the fuselage 
forward of the wings at freq. = 34.92 GHz 
(model frequency ) . 

(b) Roll plane pattern (E^) for a A/4 monopole 

above the wings 64 

28 (a) Roll plane pattern (Ee) for a KA-band 

circumferential waveguide forward 
of the wings. 

(b) Roll plane pattern (E 0 ) for a KA-band 

circumferential waveguide above the wings 65 

29 (a) Roll plane pattern (Ea) for a KA-band axial 

waveguide forward of the wings. 

(b) Roll plane pattern (E<^) for a KA-band axial 

waveguide above the wings 66 

30 The four dominant GTD terms that radiate at 

0=90°, <^=145°. 67 

31 (a) Elevation plane pattern of an axial slot 

mounted on a 4 a x 2 a prolate spheroid 
with the two-dimensional theoretical 

solution presented 69 

(b) Elevation plane pattern of an axial slot 
mounted on a 4A x 2A prolate spheroid with 
the three-dimensional theoretical 

solution presented 70 

32 Fore and aft sectors outside of which our 

"Roll Plane Model" is valid 71 

33 Illustration of data points taken from the scale 

model aircraft for the determination of best fit 
elliptic cylinder using a digital computer 74 


yi. 



Page 


Figure 


34 (a) Composite aircraft modeland coordinate 

systems illustration 76 

(b) Illustration of various coordinate systems used 
in our volumetric pattern analysis for center- 
line mounted antennas 77 

35 Illustration of belt region in which elevation 

plane model analysis is employed 78 

36 Illustration of the spherical coordinate system 81 


37 (a) Elevation plane pattern of an axial slot mounted 

on a 4Xx 2A prolate spheroid using the surface 
of revolution approach. {Geometry illustrated 
in Fig. 30) . 

(b) Elevation plane pattern of an axial slot 
mounted on a 4X x 2X prolate spheroid using 
newly developed volumetric solution. 


(Geometry illustrared in Fig. 30) 85 

38 (a) Computer simulated model for the. fuselage 

profile of a KC-135 aircraft (side view) 86 

(b) Computer simulated model for the cross^ 
section (at antenna location) of a KC-135 

aircraft (front view) 87 

(c) Computer simulated model for a KC-135 aircraft 

(top view) 88 

39 (a) Elevation plane pattern for a x/4 monopole 

mounted forward of the wings on a KC-135 
aircraft (with radome and vertical stabilizer 

included) 89 

(b) Elevation plane pattern for a x/4 monopole 
mounted above the wings on a KC-135 aircraft 
(with radome and vertical stabilizer included)... 90 


40 (a) Elevation plane pattern for a circumferential 

KA-band waveguide mounted forward of the wings 
on a KC-135 aircraft (with radome and vertical 


stabilizer included) 91 

(b) Elevation plane pattern for a circumferential 
KA-band waveguide mounted above the wings on 
a KC-135 aircraft (with radome and vertical 
stabilizer included) 92 


yii 



Figure Page 

41 (a) Elevation plane pattern for an axial KA-band 

waveguide mounted forward of the wings on a 
KC-135 aircraft (with radome and vertical 

stabilizer included) 93 

(b) Elevation plane pattern for an axial KA-band 
waveguide mounted above the wings on a KC-135 
aircraft (with radome and vertical stabilizer 
included)... 94 

42 Computer simulated model for a x/4 monopole 
mounted at station 222 on the bottom of the fuselage 
of a Boeing 737 aircraft (with radome/nose section 

being modelled by a flat plate) 96 

43 Elevation plane pattern of a x/4 monopole mounted 

at station 222 on the bottom of a Boeing 737 air- 
craft. (<f)=0° at the left; (|)=180° at the right)....... 97 

44 (a) Computer simulated model for the fuselage 

profile of a Boeing 737 aircraft (side view). 

The antenna is located at station 220 on top 

of the fuselage 99 

(b) Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 
aircraft (front view). The antenna is located 

at station 220 on top of the fuselage 100 

(c) Computer simulated model for a Boeing 737 
aircraft (top view). The antenna is located 

at station 220 on top of the fuselage 101 

45 Elevation plane pattern of a x/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

(<|) = 0° at the left; (j)=180° at the right). 102 

46 Roll plane pattern of a x/4 monopole mounted at 
station 200 on top of a Boeing 737 aircraft. 

(<f)=90° at the left; 4i=270° at the right) 103 

47 Azimuth plane pattern of a x/4 monopole mounted at 

station on top of a Boeing 737 aircraft. (e=92°) 104 

48 Illustration of the coordinate system used for 

experimental measurements... 105 

49 Elevation plane pattern of a x/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

(<))=10'’ at the left; (})= 190° at the right).. 108 


vili 



Figure Page 

50 Elevation plane pattern of a A/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

((f)=20° at the left; (})=200° at the right) 109 

51 Elevation plane pattern of a X/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

((1)=30° at the left; 4>=210° at the right) 110 

52 Elevation plane pattern of a X/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

((t)=40° at the left; (|)=220° at the right) Ill 

53 Azimuth plane pattern of a x/4 monopole mounted at 
station 220 on top of a Boeing 737 aircraft 

(9=50°) 112 

54 Azimuth plane pattern of a x/4 monopole mounted at 

station 220 on top of a Boeing 737 aircraft 
(0=60°).... 113 

55 Azimuth plane pattern of a x/4 monopole mounted at 
station 220 on top of a Boeing 737 aircraft 

(9=70°) 114 

56 Azimuth plane pattern of a x/4 monopole mounted at 
station 220 on top of a Boeing 737 aircraft 

(0=80°) 115 

57 Azimuth plane pattern of a X/4 monopole mounted at 
station 220 on top of a Boeing 737 aircraft 

(0=90°) ,. 116 

58 Azimuth plane pattern of a x/4 monopole mounted at 
station 220 on top of a Boeing 737 aircraft 

(0=100°) 117 

59 Azimuth plane pattern of a x/4 monopole mounted at 
station 220 on top of a Boeing 737 aircraft 

(0=110°) 118 

60 Azimuth plane pattern of a x/4 monopole mounted at 
station 220 on top of a Boeing 737 aircraft 

(0=120°) 119 

61 Calculated volumetric directive gain pattern of a 
x/4 monopole mounted at station 220 on top of a 

Boeing 737 aircraft 120 


i/( 



Figure Page 

62 Measured volumetric directive gain pattern of a 
A/4 monopole mounted at station 220 on top of a 
Boeing 737 aircraft. This result was performed 

at NASA (Hampton, Va.) 121 

63 (a) The scale model (1/11) of Boeing 737 aircraft 

in Anechoic Chamber at NASA (Hampton, Va.) 124 

(b) Computer simulated model for a x/4 monopole 

mounted at station 222 on the bottom of the 
fuselage of a Boeing 737 aircraft 125 

(c) Elevation pattern of monopole on a 1/llth scale 

model of Boeing 737 aircraft 126 

(d) Phase error plot of monopole on a l/llth scale 

model of Boeing 737 aircraft 127 

(e) Elevation pattern of circumferential slot on 

a 1/llth scale model of Boeing 737 aircraft.. 128 

(f) Phase error plot of circumferential slot on a 

1/llth scale model of Boeing 737 aircraft 129 

(g) Elevation pattern of a^cial slot on a 1/llth 

scale model of Boeing 737 aircraft 130 

(h) Phase error plot of axial slot on a 1/llth scale 

model of Boeing 737 aircraft 131 

64 Calculated volumetric directive gain pattern of a 
A/4 monopole mounted at station 222 on the bottom 

of a Boeing 737 aircraft (with landing gear)*--- 133 


65 Measured volumetric directive gain pattern of a 
A/4 monopole mounted at station 222 on the bottom 
of a Boeing 737 aircraft (with landing gear). 

This result was performed at NASA (Hampton, Va.) 134 

66 (a) Computer simulated model for the fuselage profile 

of a Boeing 737 aircraft (side view). The antenna 
is located at station 222 on the bottom of the 
fuselage (with landing gear and radome present).. 136 

(b) Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at station 
222 on the bottom of the fuselage (with landing 

gear and radome present) 137 

(c) Computer simulated model for a Boeing 737 aircraft 
(top view). The antenna is located at station 
222 on the bottom of the fuselage (with landing 

gear and radome present) 138 


X 



Figure 


Page 

67 Elevation plane pattern of a x/4 monopole 
mounted at station 220 on top of a Boeing 737 

aircraft (<j>=0° at the left; 4>=180° at the right) 139 

68 Elevation plane pattern of a X/4 monopole 
mounted at station 250 (off center) on top of 
a Boeing 737 aircraft (<|)=0° at the left; 

(j)=180° at the right) 140 

69 Elevation plane pattern of a x/4 monopole mounted 
at station 305 on top of a Boeing 737 aircraft 

(([)=0° at the left; 4i=180° at the right) l4l 

70 (a) Computer simulated model for the fuselage profile 

of a Boeing 737 aircraft (side view). The antennas 
are located at station 220, 250 (off-center), and 
305 on top of the fuselage, respectively 142 

(b) Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at 

station 220 on top of the fuselage 143 

(c) Computer simulated model for a Boeing 737 
aircraft (top view). The antenna is located 

at station 220 on top of the fuselage.. 144 

(d) Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at 

station 250 on top of the fuselage (off center).. 145 

(e) Computer simulated model for a Boeing 737 
aircraft (top view). The antenna is located 
at station 250 on top of the fuselage 

(off center) 146 

(f) Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at 

station 305 on top of the fuselage 147 

(g) Computer simulated model for a Boeing 737 
aircraft (top view). The antenna is located 

at station 305 on top of the fuselage 148 


xi . 



Page 


Figure 


71 (a) Computer simulated model for the fuselage 

profile of a Boeing 737 aircraft (side view). 

The antenna is located at station 250 on top 

of the fuselage (off center line) 150 

(b) Computer simulated model for the cross 
section (at antenna location) of a Boeing 737 
aircraft (front view). The antenna is located 
at station 250 on top of the fuselage (off 

center line) 151 

(c) Computer simulated model for a Boeing 737 
aircraft (top view). The antenna is located 
at station 250 on top of the fuselage (off 

center line) 152 

72 Calculated volumetric directive gain pattern of 

a A/4 monopole mounted at station 250 (off center) 

on top of a Boeing 737 aircraft 153 

73 Measured volumetric directive gain pattern of 

a A/4 monopole mounted at station 250 on top of a 
Boeing 737 aircraft (off center) with landing gear. 

This result was performed at NASA (Hampton, Va.) 154 

74 MLS Antenna space coverage diagram. (Each level 
indicates the minimum required dB level as compared 

to an isotropic source) 155 

75 Elevation plane pattern of a A/4 monopole mounted 
at station 950 on the bottom of a Boeing 737 

aircraft (<fi=0° at the left; ((>=180° at the right) 156 


76 Roll plane pattern of a A/4 monopole mounted at 
station 950 on the bottom of a Boeing 737 

airclraft. (<t>=90° at the left; (|)=270° at the right).... 157 

77 Azimuth plane pattern of a a/ 4 monopole mounted at 
station 950 on the bottom of a Boeing 737 aircraft 


(0=92°) 158 

78 (a) Computer simulated model for the fuselage 

profile of a Boeing 737 aircraft (side view). 

The antenna is located at station 950 on the 

bottom of the fuselage 159 


^tt 



Figure 

78 

79 

80 

81 

82 

83 

84 

85 

86 


(b) Computer simulated model for the cross 

section (at antenna location) of a Boeing 
737 aircraft (front view). The antenna is 
located at station 950 on bottom on the 
fuselage 

(c) Computer simulated model for a Boeing 737 
aircraft (top view). The antenna is located 

at station 950 on the bottom of the fuselage 

Calculated volumetric gain pattern of a x/4 
monopole mounted at station 950 on the bottom of 
a Boeing 737 aircraft- 

Measured volumetric directive gain pattern of a 
A/4 monopole mounted at station 950 on bottom of 
fuselage of a Boeing 737 aircraft without landing 
gear.“€-'-This Tes"ultr was performed at NASA ( Hampton, Va. ’ 

(a) Illustration of coordinate system rotation. 
(Rotate <{)o angle about z axis and then 

rotate 0 q angle about x' axis) 

(b) Illustration of coordinate system rotation 
(Rotate (})o angle about z axis and then 

rotate Qq angle about y' axis) 

Geometrical configuration of the longitudinal 
and transverse surface profile at the antenna 
location 

Front view of the roll plane (transverse) cylinder 
showing the antenna location and its surface normal... 

Geometrical configuration of the antenna location 

Illustration of the coordinate systems used in the 
roll plane model relative to the reference coordinate 
system for off center line mounted antennas 

Illustration of the coordinate system used in the 
elevation plane model relative to the reference 
coordinate system for off center line mounted antenna. 


Tciii 


Page 

160 

161 

162 

163 

168 

168 

172 

173 

174 

177 

179 


CHAPTER I 
INTRODUCTION 


Volumetric pattern analysis of fuselage-mounted airborne 
antennas at high frequencies is the object of this research. The 
primary goal of this investigation is to develop a theoretical 
solution for predicting radiation patterns of airborne antennas 
in an accurate and efficient manner. It is an analytical study 
of airborne antenna pattern problems for frequencies above about 
100 megahertz in which the antenna is mounted on the fuselage 
near the top or bottom. Since it is a study of general -type 
aircraft, the aircraft is modeled in its most basic form.’ The 
fuselage is assumed to be an infinitely long perfectly con- 
ducting elliptic cylinder in its cross-section and a compbsite 
elliptic cylinder in its elevation profile. The wing, cotkpit 
and stabilizers (horizontal and vertical) are modeled by "'n" 
sided bent or flat plates which can be arbitrarily attached 
to the fuselage. The cylinder solution employed in the analysis 
is quite adequate and useful in that reasonable approximafions 
are made in obtaining the necessary diffraction terms. This 
solution gives an accuracy of about t 1 dB in pattern compu- 
tation down to a level about 20 dB below the pattern maximum. 

Based oh the two principal plane (roll and elevation) model 
analyses from previous work [1] and extensive work on prolate 
spheroids [2], a numerical solution for the volumetric pattern 
of fuselage mounted antennas is developed. This new solution 
not only predicts the volumetric pattern accurately, it also 
reduces computer time and improves computational efficiency 
considerably over that of a previous surface of revolution 
approach [2]. The validity and capability of this new solution 
are illustrated through comparison with extensive experimental 
scale-model measurements performed at NASA (Hampton, Va.). 

One of the principal problems in the design of a reliable 
airborne antenna is the location of the antennas on the aircraft 
structure for desired radiation coverage. If modern systems are 
to function properly, the antenna pattern must meet certain 
specifications. In fact, system performance often is very much 
dependent upon the resulting patterns. 

The usual design procedure for the antenna consists of 
allocating specific locations for the antenna regardless of 
electronic system requirements but based primarily upon con- 
venience with regard to aircraft structural specifications. 

Scale model measurements, then, are used to evaluate the 
performance of the antenna system in terms of its desired 
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pattern. This approach of airborne antenna design requires a 
great deal of engineering time and money. Thus, the need for 
ah analytical investigation of the antenna system performance 
in the presence of airframe is quite apparent. 

If airborne antennas can be located on the aircraft at the 
design stage of the aircraft, then one can expect better perform- 
ance in that optimum locations and necessary structural changes 
can be anticipated. In addition, a future relocation or addition 
of antennas on an aircraft within its useful lifetime can also be 
easily accomplished through the computer simulation method. OnCe 
an optimum region is determined using computer simulated models 
of the aircraft, the antenna can be flight tested to ascertain its 
actual performance. 

Not only can these computer simulated results be used to 
determine the best location, but they can, also, determine the 
optimum antenna design for a given application. These analyses 
consider the three mutually orthogonal delta function sources 
which can be used to solve for the pattern of an arbitrary 
antenna simply by integration over the equivalent aperture 
currents [3]. 

Among the first solutions used to compute on-aircraft 
antenna patterns were the modal solutions for infinitely long 
circular [4,5] and elliptical [6] cylinders. These solutions 
modeled the fuselage by a cylinder whose elliptical cross-section 
approximated the fuselage cross-section at the antenna location. 
Arbitrary antennas were considered in these studies in which the 
antenna was mounted either on or above the fuselage. However, 
these solutions were not always adequate in that the effects of 
various scatterers such as wings and stabilizers were ignored. 

A more recent approach for solving antenna pattern and 
impedance problems is the integral equation method via moment 
methods. By enforcing the boundary conditions on aircraft 
structure, the surface currents and the resulting scattered 
fields can be found. One of the first moment solutions applied 
to aircraft problems was the wire grid technique, developed by 
Richmond [7], using a point matching scheme [8]. This solution 
required the determination of approximately 100 unknown currents 
per square wavelength in order that the wire grid adequately 
model a perfectly conducting surface. A more sophisticated 
approach has been developed by Richmond [9] in which the 
reaction technique is used to solve for the unknown currents. 

Yet this solution still requires a wire grid model of the 
aircraft with approximately 100 unknown currents per square 
wavelength. 
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Another approach is the surface-current model method [10] 
in which the surface of a conducting body is divided into patches 
with each patch having two orthogonal unknown^cpmpl ex currents. 
Using this approach, the unknown currents have been reduced to 
about 20 per square wavelength. Thus, one is able to consider 
a much larger surface using these surface patch solutions. 

However, all of these solutions are restricted to lower 
frequencies, based on the fundamental limitation on the size 
of matrices which modern computers can solve without excessive 
loss of accuracy. 

Another approach that has found great success at analyzing 
on-aircraft antennas is the Geometrical Theory of Diffraction 
(GTD). GTD is basically a high frequency solution which is 
divided into two basic problems; these being wedge diffraction 
and curved surface diffraction. The wedge diffraction solution 
has been applied to determine the radiation patterns of such 
basic antennas as parallel plate antennas [11,12,13], parallel 
plate arrays [14,15], horn antennas [16,17], parabolic reflectors 
[18,19], and rectangular waveguide antennas [20]. Both of these 
diffraction solutions have been applied in computing the patterns 
of antennas mounted on cylinders [21,22,23], rockets [24], wings 
[25,26], and aircraft [1,2,27-32]. Using this approach, one 
applies a ray optics technique to determine components of the 
field incident on the various scatterers. Components of the 
diffracted field are found using the GTD solutions in terms 
of rays which are summed with the geometrical optics terms in 
the far field. The rays from a given scatterer tend to interact 
with the other structures causing various higher order terms. 

In this way one can trace out the various possible combinations 
of rays that interact between scatterers and determine and 
include only the dominant terms. Thus, one need only be con- 
cerned with the important scattering components and neglect 
all other higher-order terms. This makes the GTD approach 
ideal for a general high frequency study of on aircraft antennas 
in that only the most basic structural features of the aircraft 
need to be modeled. i 

Since GTD is essentially a high frequency solution, the 
lower frequency limit of this solution is dictated by the spacings 
between the various scattering centers in that they should be at 
least a wavelength apart. In some cases even this requirement 
can be relaxed. Under this restriction, the low frequency limit 
is around 100 MHz. The upper frequency limit is dependent on 
how well the theoretical model simulates the important details 
of the actual structure. 


3 



The basic approach applied in the present study is based 
on previous work [1] which demonstrated the capability of the 
numerical solutions to predict the radiation patterns of fuselage 
mounted antennas in an efficient and economical way. If the 
volumetric patterns were found directly by analyzing rays on 
complex three-dimensional surfaces as done previously in Reference 
[2], the resulting numerical solution would be very complex, time- 
consuming, and uneconomical. Nevertheless, if certain assumptions 
can be made, the approach undertaken previously can be used to 
overcome these difficulties and simplify the problem a great deal. 

First, it has been shown by comparison with numerous scale 
model measurements that the roll plane model (an infinite elliptic 
cylinder with flat plate wings) can be extended to cover almost 
the complete volumetric pattern except for two conical sectors 
(fore and aft) [30]. The limitations of the roll plane model 
are due to the finite length of fuselage. Yet, the problem of 
finite length fuselage has been solved, previously, in an 
elevation plane analysis (an infinite composite elliptic cylinder). 
Furthermore, based on previous three-dimensional studies of geodesic 
rays which contribute to the radiation pattern of an antenna on 
Various prolate spheroids, one is able to combine the analysis 
of these two models to give the complete pattern. 

In addition, the cockpit/radome section and vertical 
stabilizer previously ignored are taken into account in the 
present study. The cockpit/radome section and stabilizers are 
approximated by flat or bent plates which can be arbitrarily 
attached to the fuselage. This is an improvement over the 
previous roll plane model [1]. Using this approach, the complete 
volumetric pattern can be obtained with a model consisting of a 
composite ell iptic cylinder to which are attached flat or bent 
plates. As a result of this simplified model, the solution is 
very efficient and requires little computer storage. 

The basic theoretical background on the geometrical theory 
of diffraction (GTD) is presented in Chapter II. Both wedge 
diffraction and curved surface diffraction are discussed. The 
GTD formulation for aperture radiation of antennas mounted on 
infinitely long elliptic cylinders (curved surfaces) is also 
presented in that chapter. This formulation is essential for 
fuselage mounted antennas. Chapter II also includes the near 
field scattering by a finite bent plate which is necessary for 
the determination of the scattering effects of secondary con- 
tributors such as wings, cockpit, and stabilizers. 

Chapter III describes the mathematical model of the aircraft. 
This chapter starts with a review of the two principal plane (roll 
and elevation) analyses and their application to the radiation 
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patterns. It, then, explains the way these two analyses are com- 
bined into a complete solution for the three-dimensional aircraft 
model . The computer si mu Tati 6'h method to approximate the cross 
section and elev^atipn profile of an aircraft fuselage is also 
discussed. 


Chapter IV presents the numerical results for the antenna 
patterns using these newly developed solutions. Chapter V is 
devoted entirely to airborne antenna research on the microwave 
landing system (MLS). This application is treated in detail in 
that extensive experimental scale model measurements are avail- 
able to verify the analytical solutions. -Volumetric patterns of 
both calculated and the experimental results are presented to 
demonstrate the validity and capability of this new solution. 
Finally, a summary of the present study and a discussion of 
future topics are presented in Chapter VI. 
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CHAPTER II 

THEORETICAL BACKGROUND 


As stated in the introduction, the Geometrical Theory of 
Diffraction (GTD) is ideal for a general high frequency study 
of on-aircraft antennas- This is particularly true when the 
scattering object is large in terms of the wavelength such as for 
electrically large aircraft. An aircraft shape is, in general, 
quite complex consisting of many complicated scattering structures 
To be able to obtain an accurate radiation pattern, one must take 
these various scattering structures into account. 

Based on past performance, the GTD has proven itself well 
suited to this type of analytical study. Not only does this 
approach fit nicely into a ray optics format, but it also pro- 
vides a means of analyzing the effects of three dimensional 
structures and identifying the significant contributions in the 
resulting antenna pattern. Consequenctly , the GTD is employed 
in this study to analyze the volumetric patterns of fuselage 
mounted antennas. 

The Geometrical Theory of Diffraction was introduced by 
Keller [33] as an extension of geometrical optics to include 
diffracted fields in the high frequency solution. The theory 
is based on the following postulates: 

(1) The diffracted field propagates along rays which are 
determined by a generalization of Fermat's Principle 
to include points on edges, vertices, and smooth 
surfaces in the ray trajectory. 

(2) Diffraction like reflection and transmission is a 
local phenomenon at high frequencies, i.e., it 
depends only on the nature of the boundary surface 
and the incident field in the immediate neighborhood 
of the point of diffraction. 

(3) The diffracted wave propagates along its ray so that 

(a) power is conserved in a tube (or strip of rays), 

(b) the phase delay along the ray path equals the 
product of the wave number of the medium and 
the distance s. 
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Using these postulates, one can expj;^ess the diffracted field 
in the same form,j;as^a geometrical optics’ fiel^ with solne ^ef- ^ 
ficient of prop(|rtf6nal i ty to the incident field at the point of 
diffraction. The" coefficient is determined from a canonical 
problem and is referred to as a diffraction coefficient. 

For practical purposes, the GTD can be divided into two 
categories: (1) wedge diffraction theory - to treat diffraction 

by edges and (2) creeping wave theory - to treat diffraction by 
curved surfaces. The GTD has evolved considerably from its 
original form in References [34] and [35]. Consequently, 
descriptions of the wedge diffraction problem and the diffraction 
by a curved surface are presented in the following sections. 



An asymptotic solution for the diffraction from a con- 
ducting wedge was first solved by Sommerfeld [36]. Originally, 
plane wave diffraction coefficients as presented by Keller were 
used as the sole GTD solution; however, as shown in Reference 
[37] diffraction of cylindrical waves is necessary in the treat- 
ment of antennas. 


Pauli [38] introduced the Vb function as a practical 
formulation to the solution for a finite angle conducting 
wedge. Recently, however, Hutchins and Kouyoumjian [39,40] 
presented a formulation for the diffracted field (Vb), which 
significantly improves the accuracy over that obtained from 
Pauli's form. This improved diffraction solution [39,40] 
provides superior results in the transition regions (near the 
incident and reflected shadow boundaries). It can be written 
in the form 


(1) Vg(L,<i),n) = I_^(L,(f.,n) + I^^(L,(j),n) 
where 

-j(kL + tt/4) _ /_+a\ 

(2) I^^(L,((),n) — — j ^ 

j n ^ 2 tt ' ' 

kL f°° • 2 

e"^^ ® e”^^ dx + [higher-order terms] 
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and where the higher-order terms are negligible for large kL 
and with n defined by the wedge angle [WA = (2-n)ir]. Also, 
a = 1 + cos (<j)-2n7rN) and N is a positive or negative integer 
or zero, whichever most nearly satisfies the equations 


(3a) 

2mrN-(() = -rr 

for I 

(3b) 

2nirN-(j) = +TT 

for I 


The variables L and ^ are defined later. 

The three-dimensional wedge diffraction problem is depicted 
in Fig. 1. A source whose radiated E field is given byT’(s) is 
located at point s'. It can be an arbitrary electric or magnetic 
source causing a plane, cylindirical , conical, or spherical wave 
incidence on the wedge. The diffracted vector field at observa- 
tion point s can be written in terms of a dyadic diffraction 
coefficient. Kouyoumjian and Pathak [41,42] have given a more 
rigorous basis for the GTD formulation and have shown that the 
diffracted fields may be written compactly if they are in terms 
of a ray-fixed coordinate system. The ray-fixed coordinate 
system is centered at the point of diffraction Qe (or points of 
diffraction in the case of plane wave incidence). Qe is a unique 
point (or points) for a given source and observation point. The 
incident ray diffracts as a cone of rays such that the half cone 
angle “ 6'o> the angle which the incident ray makes with the 
edge (-see Fig. 1 ) . 

The relationship between the orthogonal unit vectors associ- 
ated with the ray-fixed coordinate system (s',B' s,B.,^) are 

given by 

A 

(4) i = S'o x$' 

s — Bq X <j) , 

where I is the incident direction unit vector, and s is the 
diffraction direction unit vector. The diffracted field is, 
now, given by 

(5) E^(s) = E^Q^) • f^(s,I) A(s) 
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Where 


(6) D^(s,i) = - i'i 


Equation (6) is a dyadic diffraction coefficient as given in 
Reference [41,42]. The quantities D5 and Dh are the scalar diffrac- 
tion coefficients for the soft (Dirichlet) boundary condition at 
the surface of the wedge and the hard (Neumann) boundary condition, 
respectively. These diffraction coefficients Dg and are 
related to the Vg function described earlier by a constant as 


(7) 


Ds = 

h 


^/L 


JkL 


sin 3, 


For our purpose, it is more convenient to write the diffracted 
field in terms of the Vg function in Eq. (7) as* 


( 8 ) 




■-''b 

0 ' 


'e,; (Qj)' 



_ 0 



t 

m 


A(s) 


Sin 3 


where 


(9) Vg = Vg(L,$',n) + Vg(L,$'^,n) 


The minus sign (Vg) applies for the E-field vector parallel to the 
edge with boundary condition 

I 


*If an edge fixed coordinate system is used, Eq. (8) takes the form 
of a 3 X 3 matrix. 
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The plus sign (V 3 ) applies for the E-field vector perpendicular 
to the edge with boundary condition 


01 ) 


= 0 


wedge 


The angular relations are expressed by 


(12) - 4> + (j)' 


where the minus sign ($') is associated with the incident field 
and the plus sign ($''■) with the reflected field. The quantity 
A(s) is a ray divergence factor given by [41,42] 


(13) 


A(s) 



plane, cylindrical (s =p), 
and conical wave incidence 

spherical wave incidence 


and L, distance^ parameter, is given by [41,42] 

plane wave incidence 


(14) 


L = 


s sin^ e. 


<; 


p+ p 


* 2 

s's, sin 6 ^ 
s + s ' 


cylindrical wave incidence 

conical and spherical wave 
incidence. 


For the two-dimensional wedge problem, illustrated in 
Fig. 2, where there is cylindrical wave incidence with B =90°, 
Eq. ( 8 ) reduces to give 



e[J(p .41) 


'-''b 

0 ' 

‘E,i(p’. <(>')■ 

(15) 

_E^(p ,({))_ 


. 0 

''b-' 

_e]_(p ' , 4 >' )_ 



jk 

e 


2 ^ 

p+p' 
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In the far field (p>>p') this becomes 
(16) 

Putting this in ray form and factoring out 
(17) 


'e[|(p,(|))' 


■-V- 0 ■ 

Ej(p-,4>’) 

_E^(p,<l>)_ 


o 

<: 

.Ei_(p •,<))■)_ 


. e 


-jkp 




R‘?,(<|)) 

% 

->'b 

0 

R,‘i(<i)') 

1 

-e- 


0 

''b 

1 


£jkp 


, one obtains 


The ray form used here is given by 

p-jkp 

n / . \ c 


(18) E(p,<t.) = R(<i.) 


JT 


Thus, R(({i) is related to the far field pattern function. 

For the three-dimensional wedge problem, where there is 
spherical wave incidence, Eq. (8) reduces to give 

(19) 


'- j 

Ei|(s, Bq,<()) 


-''b 

0 

■- — 
E||(s'.s;,*') 

E 1 ( s , , (j) ) 


0 

-''a. 

Ei(s',s;,*') 


jks^s_si^ -jks 
s' e S + s' 


s+s 


In the far field (s>>s'), we have 

( 20 ) 


^1/ (s 

'X/ 

’-''b 

0 ' 

‘E,i(s*,8;,<^') 

- Ej ( s , S q j (|> ). 


- 0 

< 

.E_i(s',S’,^'] 


2 

jks 'sin 


s e 
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it is seen that 


p-jks 

Putting this in ray form and factoring out — - — 





’-''b 

0 


(21) 



. 0 


-Rl(6;,+'r 



2 

cos 3 



It is interesting to note that in the principal plane (3 q = 90°) 
the ray form of the three-dimensional case takes on the same form 
as the two-dimensional problem. 

The total ray value at the observation point s is given by 
the sum of the geometrical optics terms and the diffracted terms 

(22)< R^(s) = ■^•^•(s) + ^(s) 


where 


(23) 


R (s) + R (s) incident and reflected region I 


^. 0 . 


(s) H 



incident region II 


0 shadowed region III 


and R (s) may be determined from the image of the source term using 
basic geometrical optics techniques. These three regions are 
illustrated in Fig. 2 for a two-dimensional wedge diffraction 
problem. 


B . Diffraction by a Curved Surface 

When an incident ray strikes a smooth, curved perfectly 
conducting surface at grazing incidence, i.e., at the shadow 
boundary, a part of its energy is diffracted into the shadow 
region. To describe this phenomenon, Keller [43] introduced 
a second class of diffracted rays which is now well known as 
creeping waves. These ray paths include the points Q-j and Qp 
which form a curve on the diffracting surface, as illustrated 
in Fig. 3. However, the actual concept of creeping waves was 
introduced by Franz and Depperman [44,45]. The basic concept 
as presented in the following discussion is taken from 
"Asymptotic High-Frequency Methods" by Kouyoumjian [46]. 
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OBSERVATION POINT 



Fig. 2 —Two-dimensional wedge diffraction geometry. 
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SHADOW 

BOUNDARY 



F1g. 3— Diffraction by a smooth curved surface. 


The diffraction by a smooth curved surface is shown in 
Fig. 3 in which 0 is the source point and P is the observation 
point in the shadow region. Applying Fermat's principle, the 
line 0 QtQ 2P is the shortest distance between 0 and P which does 
not penetrate the surface. In detail, a ray incident on the 
shadow boundary at Q] divides; one part of the incident energy 
continues straight on as predicted by geometrical optics, a 
second part follows the surface s into the shadow region as a 
surface ray shedding diffracted rays tangentially as it 
propagates where f, ft, and S are the unit vectors in the 
direction of incidence, normal to the surface s and bi normal 
to the surface (6 = t x n), respectively. The incident field 


EMQl).may be resolved into 
(n • P(Q-|) and S • TCQ])) 
components induce 
pendently of each 
Q2- 


its normal and binprmal components 
It is assumed that these two 
surface ray fields which propagate inde- 
other along the geodesic arc between Q] and 
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From Reference [46] the binormal surface ray field at Qi 
is related to the bi normal component of the incident field at Q] 
by 

j ( Q 1 ) ^ 4 ' 

(24) A(Q^) e ' = D^(Q^) b^ • e'(Q^) 

where Ds(Q]) is the scalar diffraction coefficient for a soft 
surface, the amplitude of the surface ray is assumed to be 
governed by the conservation of energy between a pair of adjacent 
surface rays. Hence, the amplitude behavior of the fields is 
given as 


(25) 


A(Q2) = A(Q^) 


dn-j 



c«(t') dt' 


where 


dni and dn^ = the separation between a pair of rays at 
and Q 2 , respectively. 

a(t) = the attenuation constant wich is a function 
of t, the coordinate along the surface ray, 
because it depends on the local radius of 
curvature and its derivatives. 


The attenuation constant a(t) is introduced due to the 
tangential shedding of rays as the surface ray, propagates . It 
is seen from Fig., 3 that Q 2 is a caustic of the diffracted field 
and the second caustic is located at a distance p from Q 2 Thus, 
the binormal component of the diffracted field which radiates 
from Q 2 towards P can be found, as in the previous edge diffrac- 
tion case, with one of the caustics used as a reference point and 
is given by 


(26) 


4<t>(Q?)r 


b, . E^(P) = D,(Q,) A(«3) 


-jks 
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From Eqs. (24), (25), and (26), there results 



It is found that b-| • E^(Qv) excites an infinity of surface 
ray modes each with its own diffraction coefficient and attenua- 
tion constant. Thus, the expression in Eq. (27) is replaced by 


(28) 




g-j[k(t+s)] 



Equation (28) relates the diffracted field at P to the incident 
field at for the soft surface boundary condition. 

An expression similar to Eq. (28) is also obtained for the 
normal component of the incident field; in this case, the scalar 
diffraction coefficients and attenuation constants for the hard 
surface replace those of the soft surface. Therefore, the vector 
diffracted field at P can be written in terms of the electro- 
magnetic field incident at Q.| as 


+ bjb, 


in which v(l,2), u(l,2) are equal to ^ 

dn-| ni/f “ ‘ , “m^^ ^ 

(30) I D^(l) D„(2) e ' 

V onp . ^ m' ' m' ' 
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with the subscripts h, s, respectively, added to Dfp and 
Note that Q] and Q 2 have been replaced by 1 and 2 for the sake 
of brevity. Finding dni , dn^* and p is simply a matter of 
differential geometry involving the rays and surface; this is 
discussed at length in Levy and Keller [43]. The generalized 
diffraction coefficient and attenuation constant can be found in 
Reference [47] . 

The diffraction thus far discussed is applied to the open 
curved surface. For a closed surface, each surface ray mode 
produced at Q] encircles the surface an infinite number of times. 
The length of the surface ray path for the nth encirclement is 
t+nT where T is the circumference of the closed surface. These 
multiple-encircling rays may be summed to contribute 



-jkT - 
e 


ajf) dt’ 
m 




to the denominator of the diffracted field. It is interesting to 
note that there must be another pair of diffraction points, Q 3 and 
Q 4 , for the closed surface as shown in Fig. 4. Therefore, the 
field at any point P in the shadow region is the sum of these two 
diffracted fields from Qt~Q 2 and Q 3 -Q 4 . The total field at any 
point in the illuminated region is, by the superposition principle, 
the sum of incident, reflected and diffracted fields. A detailed 
discussion of this subject can be found in Reference [47], 

C’. Diffraction by an Infinitely 

Long Ell iptic Cyl inder 

An important special case of this GTD solution is the one in 
which the antenna is mounted directly on the curved surface. This 
is especially true for our purpose, since the radiation from slots 
and monopoles mounted on smooth curved surfaces is pertinent to 
the design of flush-mounted antennas for aircraft and spacecraft. 
Recently, Pathak and Kouyoumjian [48-51] have extended the GTD 
technique to treat the radiation from apertures or slots in convex 
perfectly conducting surfaces. This extension of GTD has been 
successfully applied to circular and elliptic cylinders, spheres, 
and spheroids. A similar GTD analysis of the radiation from 
monopoles on a convex surface has, also, been recently completed 
by Pathak and Luebbdrs [52] with the same degree of success. 
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Fig. 4—Dif fraction by a smooth closed cylindrical surface. 

The GTD solutions for analyzing the radiation of antennas 
mounted on convex surfaces are found from asymptotic solutions 
of appropriate canonical problems [48-51]. In the deep shadow 
region, the surface rays excited by the antennas account entirely 
for the field there; whereas, the geometrical optics ray field 
adequately describes the field in the illuminated region (for 
sufficiently large closed convex surfaces, the contribution from 
the surface rays is negligible). In the transition region 
adjacent to the shadow boundary, the Fock-type functions are 
employed to describe the field. This field reduces uniformly 
to the surface ray and geometrical optics solutions outside the 
transition regions. This modification in the transition regions 
is required since the ordinary ray solutions fail therein. 
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Using the GTD solution, a launching coefficient is intro- 
duced to relate the antenna field to the boundary layer surface 
waves which propagate around the surface along geodesic paths. 
Energy is continually diffracted by the surface wave in the 
tangent direction to the propagation path. This diffracted 
energy is related to the surface rays by a diffraction coef- 
ficient which is dependent on the surface geometry at the point 
of diffraction. The surface wave energy decays along the 
geodesic path in that energy is continually diffracted. This 
decay is expressed by an attenuation coefficient which is 
dependent on the surface geometry along the geodesic path. 

The GTD solutions for infinitesimal slot and monopole 
antennas mounted on an elliptic cylinder as shown in Fig. 5 
are given, with torsional effects included [49-52], by: 

A. Monopole Case 

Lit Region 


(31) E = - sin 0 ^ 0 |^F( source) 


(32) 


Transition Region 

(a) Lit side 

_ . n' X s 

E = n r 


+ b 


Z X s. 


(sgn) 


- g*(?) 


n' • s. 


(^ kPg(Q’)y ^ 


pt (Q') 




g*(c) e • F (tangent) 



X 



(a) GEOMETRY OF MONOPOLE PROBLEM 



(b) GEOMETRY OF SLOT PROBLEM 


Fig. 5— Geometry of antennas mounted on an infinitely long 
elliptic cylinder. 
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A Nl 


(b) Shadow side 


(33) 


E = 


dt^ 

pg(Q) 

dhp 

yg(Q’ 


, 1/6 


in g*(e) + 




b (sgn) -1 

g-jk**- . p (tangent) 


|l/3 


cPg(Q')yj 


9 *U) 


Deep Shadow Region 


(34) 


E = I 

3 


n. E^ + b.(sgn) 

_J 3 3 Jpt(Q') 


-1 e! 


F. (tangent) 

vJ 


B. Slot Case 
Lit Region 


(35) E = [(e-j sin g - 62 cos g) x s]- F (source) 


Transition Region 
(a) Lit side 


(36) 


1 , 

r 


n' X s„ 


1" 

sin cos g - 

0 

cos a sin g 
s 

2 



g*(5) 


+ j A(aii 

'’t(Q') 


n' • s 

0 

Ci 

2 


Vq(Q')’ 

2 

1/3 

« 

V ^ 

UPg(Q') J 


sin a„ sin 6g*(c) 


^ -J(ks^ + C^/ 3 ) 

e • F (tangent) 


y 
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(b) Shadow side 


(37) 


Y-ir'^o 

L o 


2d dtj; 


Ml 


. 1/6 


’g(Q') 


\t 'j n [sin (og -e) g*(5)] 


+ b 


Ml jr _^_ 

L^t(Q') ^ _kPg(Q') J 
• F (tangent) 


1/3 


sin ttg sin 6g*(c) 


Deep Shadow Region 


(38) 




I n. sin(a - b) + b. -3-15-) sin a. sinB E^ 
^ J s J J ‘^•"(Q') ^ ^ 


where 


F. (tangent) 


ph I nh ,h 

E “ ) D„ L e 

4 di- m=0 


dij; 1 

d^. m m 


f • (n' X z') 
sgn = ^ ‘~ 


V ■ (n' X z')| 


5 = + 




(1 

1 

1 o -? • 

for' 

J 

Q' ^9 MJ 

5 

Ik"" 


g*(‘)>g*(‘) complex conjugates of the Fock functions [48] 
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n' ,b‘ ,t' ,n,b,t normal, binormal, and tangent unit 

vectors to the surface, at source point Q 
and diffraction point Q, respectively. 


F(-) 

phase factor required to refer the phase 
of a ray to the center of the coordinate 
system 


spread factor and equals unity for this 
case [48] 

Pg>P, 

longitudinal and transverse radii of 
curvature 

2' ,T ' 

unit vectors in the principal directions 
in which the geodesics have no torsion 

^0 

unit vector pointed from the source to 
the observation point 

^0 

the distance from the source to the 
observation point 


Note that the superscripts h and s indicate the hard and soft 
boundary conditions, respectively. 

The launching coefficients are given by [48] 


"m 


m 



1/3 

A'(-q ) 
1 ' ^m' 


at the source Q' 


= 

m 




v2/3 


"I Up, 




at the source Q' 


where Dfp is defined in Table I. The subscript m refers to the 
m^h mode of the boundary layer surface wave. Thus, ym is the 
propagation constant for the mth mode surface wave such that 
Ym = am + jk» where am is defined in Table I. The quantity 
q„, is the root of Miller-type Airy function A-j(-qfn), that is, 
Ai(-qm)=0. Similarly qm is the root of A|(-qm), the derivative 
of A-j(-qni) with respect to the argument of the Airy function. 
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The incremental arc length along the geodesic path is expressed 
by di.. The summation over j in the shadow region indicates that 
several terms can contribute_in thaj^^region^ It^s noted that_for 
a slot oriented at 45° with respect to the z-axisT Eqs. (36) to 

(38) are not sufficiently accurate for calculating the radiation 
pattern as a function of ({• where e = 45°. For this special case 
and fo^r cases sufficiently close to it, an additional correction 
term should be included as indicated in Reference [50,51]. 

One must first find an efficient solution for the geodesic 
paths on the elliptic cylinder surface in order to analyze this 
problem successfully using GTD. A preferred coordinate system 
for the elliptic cylinder is illustrated in Fig. 6 and defined by 

X = d cosh u cos V = a^ cos v 

(39) y = d sinh u sin v = b^ sin v 
z = z 

where 2d is the distance between the foci of the ellipse. 

Note that for u = uf, where Uf = tanh~^ (bf/af) , the preceding 
equations define an elliptical surface. Thus, any point on 
the elliptical surface is expressed by v, which varies from 
0 to 2tt. 

Using the calculus of variations, the geodesic paths on 
an elliptical surface are given by 



Note that v^ and Vf are, respectively, the initial and final 
values of v along a given geodesic path. If one defines the 
geodesic starting direction by the angle (as) as shown in 
Fig. 5, then C = -cos as- The advantage of this geodesic 
solution lies in the fact that the integral can be quickly 
evaluated using numerical techniques. The important parameters 
of this problem are listed below: 
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Fig. 6- 


-Diagram showing the elliptic cylinder coordinate system. 



(arc length) 
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^ -a^ sin vx + cos vy 

1 I 2 2~ 2 2 

si sin V + cos v 


02 = z 

^ A yy 

t = sin a^e, - cos a 

SI s 2 


(curval inear coordinates) 


(unit tangent vector) 


b^ cos vx + a^ sin vy 

2 2 2 2 
a^ sin V + b^ cos v 


(unit normal vector) 


b = t X n = -cos a^e, - sin a 

si s 2 


(a^^ sin^ V + b^^ cos^ v)^^^ 

. 

a^b^ sin a 
f f s 


(unit bi normal vector) 


(longitudinal radius of curvature). 


Using the preceding relations, one can employ (31)-(38) to 
determine the total radiated fields in the whole space except for 
two small sectors around the cylinder axis where is near 0 or tt, 
since the solution fails in these regions. 

D. Near Field Scattering by a 
Finite Bent Plate 


The near field scattering by a finite bent plate is a 
relatively new topic at higher frequencies where the plate is 
large in terms of the wavelength. The solution presented here 
is a practical application of the three-dimensional wedge dif- 
fraction theory given earlier. The source is defined by its 
location and far-field pattern. The far-field pattern of the 
source is appropriate in that the plate is located at least 
2 d2/a away from the source where D is the maximum dimension of 
the source. The finite plate is simply specified by location 
of its n corners. The junction edge is defined by the first 
corner specified plus an additional corner (MC) defined as input 
to the computer program. The plate is initially flat (a = 180° in 
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Fig. 7). It can then be bent about the line joining corner #1 
and MC such that 90° < a < 270°. 


STATIONARY 


X SOURCE 
LOCATION 
, <x»,y5,z») 


MC+2 


MOVING 

PLATE 


y 


Fig. 7 — Bent plate geometry. 


It is known that for a given scatter direction there is only 
one point along an infinitely long straight edge at which the 
diffracted field can emanate for a near zone source. Thus, this 
point must be found for each of the n edges that describe the 
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finite plate. There are many ways of finding this diffraction 
point, one of which is described here. Since it is known that 
Bq = (sje £ig. 1) , it is obvi ous that 

(41) e„ • i = e„ . d 
'mm 

where Gm, I, and d are, respectively, the edge unit vector, 
incident direction unit vector, and diffraction direction unit 
vector. Since^the scatter direction is known (6s»<t>s)» the 
value of epi • d = c^ is easily computed at each edge. Ope needs 
only search along the edge to find the point where e^ • I = c^. 

Once the diffraction point is located, one must find the 
diffracted field value from the mth edge. The far field pattern 
of the source can be written as 

(42) E^(e,(j)) = [0 F(9,({)) + <{)G( 0 ,i))] r — = R(0,i>) i — 


where s' is the range from the source to the field point. Using 
the geometry illustrated in Fig. 7 and applying the results 
presented earlier, one finds that 


(43) 



-Vr 


K 


-j[k(s' - y) -kpp] 


L 0 


'B-> 


L- J 


where 




R(e,.+,) • 

R(9^-J<t>.j) ' <t>Q 
2 

ks sin 


sin 0 ^ cos <t>^ + y. Sin 0^ sin it>^ + z.„ cos 0^ 

dp s s "^dp s s dp s 


Vg(kp.p/k^, 2) + Vg(kp.p/k , <|)+<|)', 2)^ 


I x,3^' 
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The coordinates (x^p, y^p, z^p) define the point of diffraction. 

Once these terms are determined, the total diffracted field 
in ray form from a general m^h edge is given by 

(44) * 3 ) ■ 


where = d x Bq- Using the superposition principle, the total 
singly diffracted field in ray form by the n edges of the plate 
is given using Eq. (44) by 

(45) R“(e^, *^) = I * 3 ) . 


The first step in calculating the reflected field is to find 
the locations of the image sources, which are uniquely determined 
once the planes of the flat plates are defined relative to the 
source location. In fact, the image is located along lines which 
are orthogonal to the different portions of the plate and posi- 
tioned an equal distance on the opposite side of the plate. 

With the image position known, one needs to determine if the 
reflected field contributes to the total scattered field using 
the geometrical optics approach. If the reflected field is a 
contributor, the ray from the image source in the scatter direc- 
tion (85, <J>s) must pass through the finite plate limits. Thus, 
one must find the location of the intersection point of this ray 
and the plane containing the flat plate. This can easily be 
accomplished using vector analysis. One can, then, predict 
within certain limits whether this intersection point falls 
within the bounds of the finite flat plate. 

If reflections do occur, the reflected field from the image 
source can be written in ray form as 

(46) ’^’s^ 

jk[x^. sin cos + y. sin sin <()g + z^. cos e^] 


where e and <() are related to the image source coordinate system 
with the image location defined by (x-j, y^ , z-j)- The functions 
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[F^( 0 o, 4 >s) and G’^(e 5 , (|)c)] are found by employing the boundary 
conditions on the perfectly conducting flat plate. The total 
scattered field from the_flat .plate_ls.,_then.,-.gi.v.en_by___„_,.^__ 

(47) R^( 0 ^, < 1 ,^) + r'^( 0 ^, (J,^) , . 

The four basic terms included in the present solution are as 
follows: 

(a) single diffraction of incident field as shown in 
Fig. 8 (a) 

(b) single reflection of incident field as shown in 
Fig. 8 (b) 

(c) double reflection of incident field as shown in 
Fig. 8 (c) 

(d) single diffraction of reflected field as shown in 
Fig. 8 (d). 

Each of the terms has been illustrated in a two-dimensional view 
just for simplicity in illustrating the mechanisms;, whereas, the 
actual solution is for the three-dimensional geometry. These 
terms have all been incorporated in a general bent plate solu- 
tion. Note that only those terms are included which are not 
shadowed by another portion of the bent plate. 

The accuracy of this solution is illustrated by the example 
shown in Fig. 9 where it compares very favorably with measured 
results taken for a A/2 dipole illuminating a flat plate. This 
solution is also compared with results obtained using a moment 
method patch technique [50] as illustrated in Figs. 10 and 11. 
Note that for these comparisons the plate dimensions are quite 
small in terms of the wavelength, which accounts for the small 
discrepancies in the patterns between the two solutions. This 
structure will be incorporated into the aircraft model of 
Chapter III in order to account for bent wings , moving flaps, 
vertical stabilizer, etc. 
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Fig. 8a— Single diffraction of incident field. 
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Fig. 8b— Single reflection of incident field. 
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Fig. 8c— Double reflection of incident field. 


SOURCE 

LOCATION 



33 







Fig. 10b— E-plane pattern. 
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Fig. lib— E-plane pattern. 
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CHAPTER'! I I 


MATHEMATICAL MODELING OF AIRCRAFT 
FOR PATTERN COMPUTATIONS 


In this chapter, the mathematical modeling of an antenna 
on an aircraft is described. This mathematical model is intended 
to be as simple as possible so that it can be used to simulate a 
wide variety of aircraft structures, yet accurate enough that the 
computed radiation patterns are comparable to measured results. 

The model is systematically increased in complexity to resolve 
discrepancies between calculated and measured results, yet basic 
simplicity is retained. One of the restrictions in the model is 
that the antenna under consideration is located near the top or 
bottom of the fuselage surface. 

This chapter begins with a brief review of the analyses 
using the elevation plane model and the roll plane model. This 
is followed by solutions for the volumetric patterns of antennas 
on a three-dimensional aircraft structure. The approach used to 
model the aircraft by two elliptic cylinders is explained in detail. 
The numerical technique to generate the necessary elliptic cylinders 
is also discussed. Finally, the radiation patterns in principal 
planes for various cases are calculated to illustrate the versa- 
tility of the newly developed solutions. 


A. Analysis of Elevation Plane Model 


In most cases, the dominant structural effect in the 
elevation plane is the profile of the aircraft for fuselage 
mounted antennas. An aircraft fuselage is usually a convex 
body that cannot be completely described by simple mathematical 
equations. Hence, in practice, an aircraft shape is often 
specified by points. Consequently, a method called "section 
matching method" has been developed [ 28 ] to handle the problem of 
an antenna mounted on a fuselage surface of general shapie. 

However, this method is restricted to principal plane pattern 
calculations. For off-principal plane pattern calculations 
the method is hard to apply, since a complete three-dimensional 
description of the fuselage surface is difficult to obtain. 
Furthermore, the lack of information on the geodesic paths and 
torsion effect of the curved surface create serious problems in 
the prediction of antenna patterns. Consequently, a simple 
model which simulates a general fuselage profile has been adopted. 
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Since an aircraft fuselage is usually long and slender, a 
composite elliptic cylinder appears to be a good model. This 
model consists of two semi-elliptic cylinders mounted back-to- 
back as shown in Fig. 12. Note that as discussed in the last 
chapter, the model employed here is not limited to the determi- 
nation of the elevation plane pattern only. It can be used to 
determine off-elevation plane radiation patterns also. 

Recall that the solutions presented in the last chapter are 
based on elliptic cylinder models; however, one of the advantages 
of GTD is that it can be extended to new structures after making 
certain assumptions. In this case, it is assumed that there are no 
diffractions from the junction lines of the two ellipses. This 
assumption is justified since these junctions are non-existent in 
the actual aircraft profile. Note that the GTD solution in the lit 
region does not depend on the surface parameters in that it is 
assumed the source is mounted on an infinite ground plane tangent 
to the surface at the source point. On the other hand, the transition 
and deep-shadow region solutions are modified due to their dependence 
upon the surface parameters. This modification requires that one 
use bf = b' for rays traveling to the right of the junction and 
bf = b for rays traveling to the left of the junction. The parameters 
b and b' are illustrated in Fig. 12. 



FIRST ANTENNA LOCATION 


* y 



Fig. 12 — Theoretical model of KC-135 aircraft. 
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Since the model employed is a composite elliptic cylinder, 
the surface curvature at the junction line is not continuous due 
to the different dimensions of the elliptic cylinders^ To ensure 
the assumption made earlier that no diffraction occurs from the 
junction line, the solution presented in Section (C) in Chapter 
II needs to be modified to correct for the difference in curva- 
ture when the surface wave crosses the junction line. Since no 
GTD solution for this problem is available at present, an alternate 
way to account for the continuity of the surface wave crossing the 
junction line is necessary. It is noted that due to the large radius 
of curvature at the junction line in our model the reflection 
effect is so weak that it can be neglected. A way to treat the 
continuity of the surface wave across the junction is investigated 
next. 


Consider a source located on a composite elliptic cylinder 
as shown in Fig. 13. The field diffracted by a curved surface 
can be written, according to Eq. (34) for hard boundary con- 
ditions and with only the firsJt mode retained, as 


(48) - E = CL(Q') D(Q) e 


(a + jk) ds 
J n ' 


where C is a complex constant. This result indicates that the field 
diffracted from the curved surface is related only to the launching 
coefficient at the source point and the diffraction coefficient at 
the diffraction point. Thus, when the diffraction point is on the 
right half elliptic cylinder, the field is given by 

- I ^ (a + jk) ds 

(49) E = C L(a,b) D 2 (a,b') e ; 


note that a and b are the semi minor and major axes of the left 
half ellipse and a and b' are the semi minor and major axes of 
the right half ellipse, respectively. Similarly, the field 
diffracted from point Q on the left half ellipse is expressed as 

^ (a + jk) ds 

(50) E = C L(a,b) D^(a,b) e . 


41 



Q » 

UJ 



42 


Fig. 13— Illustration of surface waves propagating across the junction line 
of a composite elliptic cylinder. 


At the junction, Eqs. (49) and (50) become 


Ej = C L(a,b) D 2 j(a,b' ) e 


(a + jk) ds 


J ^ (a + jk) ds 

Ej = C L(a,b) D^j(a,b) e 


respectively i where the subscript J denotes the junction. 

It is noted that Eqs. (51) and (52) do not agree^with each 
other which indicates a discontinuity exists in the diffracted 
field. However, this is not true in practice, since the field 
must be continuous as the diffraction point crosses the junction. 
Comparing Eqs. (51) and (52), one notes that if a modification 
factor 

rD,,,(a,b) 1 


is introduced, the field is continuous as the diffraction point 
crosses the junction. Thus, Eq. (48) is modified to give 


(53) E = C 


Dij(a,b) 
D2j(a,b' ; 


L(a,b) D 2 (a,b') e 


(a + jk) ds 


Notice that Eq. (53) is valid only for the surface wave which 
crosses the junction, i.e., the diffraction point and source 
point occur on different elliptic cylinders. 

From Table I, one sees that the factor introduced in Eq. 
(53) can be simplified to give 


Dlj(a,b) Pgj(a,b) 

D2j(a,b' ) ■ Pgj(a,b' i 


which depends only upon the ratio of the two semi>-major axes 
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A1(-x) ■ * I coi( ^ - Kt) dt. 

At‘(-K) denotes the derivative of A1(-s) with respect to the argunent 
of the Airy function. 

















of the composite ellipse. For a continuous elliptic cylinder, 
this factor reduces to unity as expected. Note, that, Eq. (53) 
can be generalized to describe the field diffracted from an 
arbitrary point on the curved surface as 

rQ 

(a + jk) ds 


(55) 


E = C L(Q') D(Q) p e 


Q' 


where 

(56) 


1 


P 


if Q and Q' are on same elliptic cylinder 

1/6 


(ell ipse with source) 

Pqj ( ell ipse without source) 


otherwise. 


In a similar manner, this factor can be determined through 
the field expression in the transition region as presented in 
Chapter II and is found to be the same as Eq. (56). 

To illustrate the newly modified solution, the elevation 
plane radiation patterns for antennas mounted on the KC-1 35 air- 
craft are now analyzed. Since the antennas of interest are located 
on top of the fuselage and along the center line, the most significant 
effects on the pattern result from the upper surface of the 
fuselage profile. The structure used to simulate the l/ 25 th 
scale model consists of a 3.75" by 60" right semi -elliptic cylinder 
and an 3.75" by 8" left semi -el 1 iptic cylinder as shown in Fig. 12. 
Using Eqs. (31) to (38) in Chapter II with the modifications just 
described, the radiation patterns in the elevation plane for 
antennas mounted on our model can be easily determined. 


The elevation plane patterns for a A/4 monopole mounted 
forward and over the wings are illustrated in Fig. 14 (a) and (b), 
respectively. The patterns for a circumferential KA-band wave- 
guide are illustrated in Figs. 15 (a) and (b). Finally, the 
results for an axial KA-band waveguide are shown in Fig. 15 (a) 
and (b). The KA-band waveguide aperture fields are simulated in 
our model by an array of infinitesimal elements as shown in Fig. 
17. 
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Fig. 


14a — Elevation plane pattern for a A/4 monopole mounted 
forward of the wings on a KC-135 aircraft. (No 
radome and vertical stabilizer included). 



180 ® 


Fig. 14b — Elevation plane pattern for a a/ 4 monopole mounted above 
the wings on a KC-135 aircraft. (No radrane and vertical 
stabilizer included). 
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180® 

Fig. 15a~Elevation plane pattern for a circumferential KA-band 
waveguide mounted forward of the wings on a KC-135 
aircraft, (No radome and vertical stabilizer included). 



180® 

Fig. 15b— Elevation plane pattern for a circumferential KA-band 
waveguide mounted above the wings on a KC-135 aircraft. 
(No radome and vertical stabilizer included). 
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180 ® 

Fig, 16a--Elevation plane pattern for an axial KA-band wave- 
guide mounted forward of the wings on a KC-135 air- 
craft. (No radbme and vertical stabilizer included) 



180 ® 


Fig, 16b — Elevation plane pattern for an axial KA-band wave- 
guide mounted above the wings on a KC-135 aircraft. 
(No radome and vertical stabilizer included). 
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- 0 . 112 “ - 0 . 056 " 0 “ 0 . 056 " 0 . 112 " - 0.0467 0 “ 0 . 0467 * 


Fig. 17— Approximation of an open-end waveguide by an array 
of 5 infinitesimal broad wall elements with dif- 
ferent weights and an array of 3 infinitesimal 
narrow wall elements with uniform weight. 
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The comparison between calculated and measured results is 
encouraging, especially when the antennas are located above 
the wings. The discrepancy displayed in the fore and aft regions, 
particularly for the case of antennas located forward of the wing, 
is due to the effects of the cockpit and nose section and the 
vertical stabilizer which are neglected in the present theoretical 
model . 

In the preceding analysis, the effect of the aircraft in the 
nose of the radome was neglected. In order to simulate a fuselage 
with a radome, a truncated composite elliptic cylinder model as 
shown in Fig. 18 is adopted. This is based on experimental data 
that the radome effect appears to be a diffraction from the bulk 
head as if the radome were not present. 

Before preceding to solve this problem, by careful consider- 
ation it can be converted into a simpler problem. With the fact 
that a slot or monopole radiating in the presence of a composite 
elliptic cylinder can be considered as an antenna itself, the 
present problem can be approximated by the radiation problem of 
an equivalent antenna mounted on a wedge type structure as shown 
in Fig. 19. where the antenna is mounted on the bottom of the 
fuselage. This equivalent antenna radiates the same pattern as 
that of a slot or monopole mounted on our previous model without 
the radome considered. The wedges are formed and defined by the 
tangent planes at the discontinuities and Q 2 . Thus, the 
radiation problem of a fuselage model with a radome included 
is transformed to a wedge diffraction problem, and can be solved 
using standard GTD techniques. 


Figure 19 illustrates the necessary components of the total 
electric field in the different regions given by 
(• 

iT-(j) ,<0<ir, (j)=0 

eq 0)1 ^ 

TT-(b , <0<Ti,d)=Tr 
0)1 

Note that Epg is ;the radiated field of the equivalent antenna 
descjribed above, E^ is the singly diffracted field from edge 1 
and E°° is the diffracted field from edge 2. It is noted that 
the contributions from any other higher order diffraction terms 
are so small compared to single and double diffraction that 
they are ignored for our purposes. 
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Fig. 18— Truncated composite ellipse as simulation of fuselage 
wi th radome i ncl uded . 






E +• E* 
eq ^ ^ 


Fig. 19— Total field distribution. 
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Using Eq. (57), the elevation plane patterns for the space 
shuttle as shown in Fig. 20 with the radome included are calcu- 
lated. Figures 21 to 23 illustrate the radiation patterns for both 
vertical (monopole and circumferential waveguide) and horizontal 
(axial waveguide) polarization sources mounted on the fuselage of 
the 1/35^*^ scale model of the space shuttle with the radome considered 
The structure used to simulate the l/35tfi scale model consists 
of a 2.55" by 12.5" left semi -ell ipse and a 2,55" by 60.0" right 
semi-ellipse. The size of the radome is 6.5" and the source is 
mounted 2" away from the radome. Figure 24 shows the radiation 
pattern of the space shuttle with radome size being 2" and 
antenna mounted 1/2" away from radome. Note that the waveguide 
used here is the same KA-band waveguide used previously. The 
experimental results, which were taken at NASA (Hampton, Va.) at 
a frequency of 35 GHz, are also presented. The agreement between 
calculated and measured results illustrates the applicability of 
the numerical analysis technique being developed in this study. 

Note also that although the elevation plane model can be used 
to compute the elevation plane pattern as well as off-elevation 
plane patterns, it is found that for our fuselage model, this 
model can only be used to compute the off-elevation plane patterns 
accurately up to approximately t30° from the elevation plane. 


B. Analysis of Roll Plane Model 

The basic aircraft to be analyzed in this section is composed 
of flat plates attached to an infinitely long elliptic cylinder. 
Originally, a circular cylinder was used to represent the fuselage 
[27], and the modal solutions were employed to determine the 
field. However, an aircraft cross-section is not circular in 
general. To be able to predict the radiation pattern more 
accurately, one has to consider a better model for the fuselage 
approximation. 

Since the roll plane cuts orthogonally across the fuselage, 
one should expect the fuselage cross-section to have a strong 
effect on the roll plane pattern. On the other hand, an aircraft 
fuselage is normally long and slender, such that its finite length 
effects are generally secondary. Consequently, the infinitely long 
elliptic cylinder representation of the fuselage for roll plane 
calculations appears to be a reasonable approximation in most 
cases. Since the antenna can be arbitrarily positioned on the 
fuselage with respect to the wings, one must consider the width 
of the wing as well as its length in order to obtain a practical 
analytic model. In order to accomplish this, the near field 
bent plate scattering solution is adapted to this new model 
such as illustrated in Fig. 25. Note that the wings are 
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Fig. 20b— Theoretical model of l/35th scale model of space 
shuttle. 
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Fig. 21— Elevation pattern of monopole on a l/35th scale model 
of space shuttle (with radome included). 
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Fig. 22— Elevation pattern of circunrferential slot on a 
l/35th scale model of space shuttle (with radome 
included). 


56 


TAIL 



NOSE 


HORIZONTAL POL. SPACE SHUTTLE 

WAVEGUIDE WITH RADOME 


(8/2") (6/2”) 

0=o» 



180 ® 

ELEVATION PATTERN 


Fig. 23— Elevation pattern of axial slot on a 1 /35th scale 
model of space shuttle (with radome included). 
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Fig. 24— Elevation pattern of axial slot on a l/35th scale 
model of space shuttle (with radome being 2"). 
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TAIL 


Fig. 25a— Reflection problem in x-y plane. 



Fig. 25b— Diffraction problem in y-z plane. 




assumed flat here and each wing can be located arbitrarily 
with any number of edges. The wings, also, can be mounted up 
or down from the central location provided that the wings are 
horizontal. 

The model, now, consists of an infinitely long elliptic 
cylinder fuselage to which finite flat wings are attached. The 
various configurations analyzed are shown in Fig. 26 looking 
from the front of the aircraft with the antenna mounted in 
each case above the wings for the models illustrated. Using 
these models, one should be able to analyze a wide variety of 
aircraft shapes. This is verified by a comparison of results 
taken on actual aircraft scale models and presented later. 

Let us first find the effective source location for the 
reflected field. Recall that in the flat plate result, the 
source was imaged and the reflected field added to the total 
solution provided the image ray passed through the finite 
flat plate (wing) limits. One must initially determine the 
effective source position and then the reflected field. With 
the source mounted on an infinitely long elliptic cylinder, 
the surface rays from the source propagate around the cylinder 
along geodesic paths, from which energy is continually diffracted 
tangentially. Now let us assume that the source does not illumin- 
ate the right wing directly (as illustrated in Fig. 25(a)) and 
proceed to determine the unique geodesic path that diffracts 
energy from a known tangent point which is then reflected off 
the wing in the desired radiation (or scatter) direction. 


The effective source position for reflections from the right 
wing in terms of the radiation direction (e^, <j>g) is given by 

= a^ cos v„ 
e f e 

(58) y^ = b^ sin Vg 



where Ve = tan” (bf/af cot (f>s). These coordinates can, then, 
be used in the flat plate problem as the effective source loca- 
tion. Note that as the desired radiation direction is varied, 
the effective source location changes. In addition, if the 
source directly illuminates the wing for a given reflection 


term, then the effective source location is simply the actual 
source location (uf, Vg-, Zgg) • A result similar to Eq. (58) 
can be found for the reflections from the left wing. Finally, 
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Fig. 26— Fuselage and wing geometry for theoretical aircraft 
model looking from the front. The antenna is always 
mounted on top of the models. 


the actual source field value used to compute the reflected 
term is determined from the 6TD solutions in Chapter II. 

Using a similar technique, the effective source locations 
for the diffracted field components may be found. The bent 
plate solution uses a search technique to find the diffraction 
point by computing the diffraction angles at selected test 
points along a given edge. Once a test point (x<j, yj, zj) is 
specified along the edge, one can find the effective source 
location (xe, ye> Ze) rising the geometry illustrated in Fig. 
25(b). Again it is assumed that the source does not directly 
illuminate the test point. One finds that the effective source 
is given by 


X. = 


2 2 

afbfXd 




2 2^ 


2 

^f^d 


afbf 


2 2 

"f^d ^ 
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(59) 


^e = 


2 2 
3fbfyd 


VdH 


2 2 

^f^d + 


k 2 2 

*^f^d 


2 2 2 2 

44 * 44 


2 . 2 

3(,bf 



Ve'v 


where 




® I 2~ 2 2 2 

J a^ sin V + b^ cos v dv. 


so 


r~2 2 2 2 

^v ~ J^f '^e '^e* 




Given the effective source location for the chosen test 
point, the search technique is applied to find the actual dif- 
fraction point along a given edge. Note that once the actual 
diffraction point is determined, the effective source of the 
diffracted field is specified by Eq. (59), and the source field 
value is, again, computed using the 6TD solutions. 
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The total field is found by summing the directly radiated 
field with the scattered field using the , superposition^ principle. 
The results for a x/4 monopole on the fuselage of a l/25th scale 
model of a KC-135 aircraft, as shown in Fig. 12, forward and over 
the wings are shown in Figs. 27 (a) and (b), respectively. The 
results for a KA-band circumferential waveguide forward and over 
the wings are shown in Fig. 28. The results for a KA-band axial 
waveguide forward and above the wings are shown in Fig. 29. The 
fuselage, in this case, is approximated by a 3" x 3" infinitely 
long elliptic cylinder. The waveguide antenna is modeled as in 
the previous section, and the agreement in each case is very 
encouraging. Again, it should be stressed that the roll plane 
model solution not only can predict. the radiation pattern 
accurately in the roll plane, it also can be extended to cover 
almost the complete volumetric pattern as shown in Ref. [30]. 


C. Three-Dimensional Model 

Approximation of an Aircraft 

The analytical solutions described in the previous sections 
provide a useful, efficient, and economical way for the evalu- 
ation, location, and design of fuselage mounted antennas based 
on their pattern performance in principal planes. However, if 
modern systems are to function properly, the antenna pattern 
must meet certain specifications. These specifications are 
usually given in terms of a coverage diagram for a particular 
sector in space which is to be met with antenna mounted on the 
aircraft not on a finite ground plane where the original antenna 
design was made. Thus, the desire for an accurate solution for the 
complete pattern performance of antennas mounted on a complex air- 
craft structure for given applications requires a more thorough study 
of ways to handle the volumetric pattern. 

If this problem is attacked directly by analyzing rays 
on complex three-dimensional surfaces as done previously in 
Reference [2], the resulting numerical solution would be very 
complex, time-consuming, and uneconomical. Nevertheless, if 
certain assumptions can be made, the approach undertaken in 
the previous sections can be used to overcome these difficulties 
and simplify the problem a great deal. 

First, it has been shown by numerous scale model measure- 
ments that the roll plane model can be extended to almost 
cover the complete volumetric pattern except for two conical 
sectors (fore and aft) [30]. The limitations of the roll plane 
model result are due to the finite length fuselage. Yet, the 
finite length fuselage has been solved, previously, in the elevation 
plane model analyses. Furthermore, based on the previous three- 
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Ftg, 27a— Roll plane pattern (E.) for a 1/25 scale model of a 
KC-135 with a A/4 monopole on the fuselage forward 
of the wings at freq, = 34.92 GHz (model frequency). 



180® 

Fig. 27b— Roll plane pattern (E ) for a x/4 monopole above 
the wings. 9 
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180“ 


Fig. 28a— Roll plane pattern (E.) for a. KA-band circumferential 
waveguide forward of the wings. 



180“ 


Fig. 28br-Roll plane pattern (Eg) for a KA-band circumferential 
waveguide above the wings. 


65 







F1g. 29a— Ron plane pattern (E.) for a KA-band axial waveguide 
forward of the wings. , 



ISO® 


Fig. 29b— Roll plane pattern (E^) for a KA-band axial waveguide 
above the wings. 
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dimensional studies of geodesic rays which contribute to the 
pattern of an antenna on various prolate spheroids, one is 
able to combine the analyser of these two pfihc^'par pTane^mb'deTs 
to give the complete pattern. In this study, computer simulation 
models were considered that would resemble a wide variety of 
aircraft shapes and yet could, also, be analyzed with reason- 
able accuracy and economy. In this case, it is quite obvious 
that the three-dimensional nature of the fuselage must be 
modeled if one is to adequately determine volumetric patterns. 

This resulted in the development of a general surface of 
revolution model for the fuselage as presented in Reference 
[2]. Through an extensive study of geodesic paths on a surface 
of revolution, the number of dominant rays that contributed to 
the radiation pattern were shown to be finite except for more 
spherical shapes. Furthermore, the computer result showed that, 
for a prolate spheroid, the dominant rays needed to be con- 
sidered would not exceed four rays; in most cases, it is even 
less than that. These four rays are illustrated in Fig. 30 
in which two rays are propagating along the cross-section of 
prolate spheroid; the other two are propagating along the profile. 



Fig. 30— The four dominant GTD terms that radiate at 
e = 90°, <t> = 145°). 


67 



To demonstrate the significance of the fOur-ray contribution, 
the elevation plane pattern of an axial slot mounted on a prolate 
spheroid was calculated using a two-dimensional (two rays) and 
three-dimensional (four rays) solution as shown in Fig. 31. 
Experimental results are also shown to verify the calculated 
three-dimensional result. It is immediately obvious that the 
back lobe region is not calculated with sufficient accuracy 
using the two-dimensional result. However, the three-dimensional 
solution is in good agreement with the measured pattern. This 
leads to a new approach to handle the volumetric patterns for 
fuselage mounted airborne antennas in a simplified and economical 
manner. 

As determined previously, the roll plane model (infinite 
elliptic cylinder with flat plate wings) can be extended to cover 
most of the volumetric pattern except for the fore and aft sectors 
(about 15° conical sector as shown in Fig. 32). To cover these 
sectors, one must incorporate the elevation plane model analysis 
with some modifications such that the effects of wings and stabilizers 
are also included. In doing so, two different elliptic cylinders are 
required to analyze the volumetric patterns; one being the cross- 
section (roll) cylinder and the other being the profile (elevation) 
cylinder. This requirement is necessary since our solutions are 
based on an infinitely long elliptic cylinder, in which the z-axis 
coincides with the axis of the cylinder. 

As discussed in the previous chapter, the most significant 
effects oh the radiation pattern result from the surface 
geometry nearest to the antenna. The curvature of the surface 
in the vicinity of the antenna location plays a dominant role 
in predicting the radiation pattern. Thus, the elliptic 
cylinders used to represent the fuselage profile and cross- 
section need to model the aircraft structure as accurately 
as possible near the antenna location. Once these two 
elliptic cylinders are obtained, one is able to proceed to 
solve for the complete volumetric pattern. 

To determine the elliptic cylinders necessary to simulate 
the fuselage profile and cross-section, one has to find some 
way to approximate the curved surface near the antenna location. 

One way to achieve this goal is to find an analytical function 
which can best approximate the known curved surface profile 
(either in elevation or roll plane) in a least mean square error 
sense. This leads to the development of a best fit routine 
to generate the necessary ellipses for the fuselage model 
through the use of a digital computer. 
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Fig. 31a— Elevation plane pattern of an axial slot mounted 
on a 4A X 2 \ prolate spheroid with the two- 
dimensional theoretical solution presented. 
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Fig. 31b— Elevation plane pattern of an axial slot mounted on 
a 4A X 2A prolate spheroid with the three-dimensional 
theoretical solution presented. 
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The theory behind the best fit routine is that a function 
is to be found to best approximate a set of points in a least 
mean square error sense. For the aircraft model, a best fit 
ellipse is desired. The mathematical expression for an ellipse 
is 

( 60 ) 4 ^ 4 ='' 

for its origin located at (x = 0, y = 0). The parameters a and 
b are semi major and minor axes of an ellipse, respectively. To 
simplify the mathematical expression, Eq. (60) can be written as 


(61) AX + BY = 1 

where X = x^,Y=y^,A = -^, and B = ^ . Let (x^. , y^. ), 

a b 

i =1, ••• n, be n points on which a best fit ellipse is to be 
generated. Substituting these points into Eq. (61), one obtains 
a set of n linear equations. 

AX^ + BY^ = 1 

AX2 + BY2 = 1 


AX„ + BY„ = 1 
n n 

In matrix form, these n equations become 


(62) ZC = I 

where 
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arid 


1 

1 


'V 

By multiplying both sides of Eq. (62) by Z, the transposed 
of Z, one obtains a simple 2x2 matrix as given by 


(63) 


or 


TIC = 1 I 


n 

o 

n 
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n 
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1 

m=l 
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1 

m=l 

m m 
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m=l 
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X 

n 
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I 

m=l 

m m 

I 

m=l 

m 


B 


I 

m=l 

_ 


which can be simply solved for A and B in the least mean, 
square serise [56], Thus, the ellipse parameters are defined. 

Before the necessary ellipse or composite ellipse can be 
determined through the numerical process, the data points that 
described the surface of the profile or cross-section oflthe 
actual aircraft fuselage must be generated. To do so, aireference 
Cartesian coordinates system is needed on a scale model drawirig 
of the aircraft. These coordinates can be best located by align- 
ing one of the axes with the center line of aircraft fuselage 
with the origin being arbitrarily chosen according to convenience. 
After the coordinate system is fixed, the positions of data points 
on the profile or cross-section can be measured from the scale 
model relative to the reference origin. The data points are 
taken in such a way that more points are needed around the 
antenna location and less points away from the source. This is 
due to the fact that the surface profile is dominant near the 
antenna locationi as described earlier. Figure 33 illustrates 
the way data points are taken from a fuselage profile. By feeding 
these data points into the best fit routine and adjusting the 
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Fig. 33--niustration of data points taken from the scale 
model aircraft for the determination of best fit 
elliptic cylinder using a digital computer. 
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origin of the coordinates in the routine, an ellipse is found to 
best fit these data points. In the same manner, the ellipse to 
approximate the cross-section of aircraft fuselage at-the antenna 
location can be obtained. Thus, the elliptic cylinders necessary 
to simulate the fuselage profile and cross-section are determined. 

Now let us consider how these two finite elliptic cylinders 
can be used to describe the volumetric pattern for fuselage mounted 
airborne antennas. First, a referente coordinate system is needed 
so that the solutions based on our two-dimensional analysis in two 
different coordinate systems can be incorporated into a complete 
solution. Since an aircraft fuselage is usually long and slender, 
its finite length effects are limited to the small sectors off the 
nose and tail. Consequently, the two rays which play the most 
significant role normally come from the cross-sectional elliptic 
cylinder. Thus, it seems more natural to have the reference 
coordinate system correspond to the roll plane model coordinate 
system, i.e., the Z^gf ^ixis is pointed aft and the X^gf axis is 
pointed vertically upward as shown in Fig. 34. Notice that the 
origin of the composite ellipse in the elevation plane is chosen 
as the origin of the reference coordinate system for convenience. 

Recall that the roll plane model is not valid in the two 

15° conical sectors in the fore and aft directions. In order to 

overcome this handicap, a belt region around the antenna source 
is chosen such that in this region the solution is obtained using 
elevation plane model analysis. Everywhere outside the belt region, 
the roll plane model analysis is used to analyze the radiation pat- 
tern as done previously. Figure 35 shows the regions in which the 
elevation and roll plane cylinder solutions are used. The angle 
a is chosen such that the roll and elevation cylinder solutions 
blend smoothly together. In fact, a is a function of the size of 
elevation and roll plane cylinders used to simulate the aircraft 
fuselage. For most cases in our model, the angle a is set at 

20°. This 20° belt has been tested and found to be satisfactory, 

based on comparisons with measured results, as will be shown later. 

As seen in Fig. 34, the elevation plane model coordinate sys- 
tem is such that the Zgigy axis points in the Y^gf direction. In 
other words, the elevation plane model coordinate system is simply 
obtained by rotating the XYZy^gf coordinate system by a 90° angle 
about the X^ef axis. Notice that the antenna here is assumed to be 
mounted on the center line of the fuselage for simplicity, but this 
is not a basic limitation. In order to determine the electromagnetic 
fields in a given radiation direction (0o ref> ')>n ref) within the 
belt region, the solution obtained through the elevation plane 
model analysis is expressed in terms of the reference coordinate 
system such that it can be superimposed with the roll plane 
cylinder solution. This requires a coordinate transformation 
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Fig. 34a--Composite aircraft model and coordinate systems illustration. 
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Fig. 34b— mustration of various coordinate systems used in our volumetric 
pattern analysis for center-line mounted antennas. 



Fig. 35— niustration of belt region in which 
elevation plane model analysis is 
employed. 
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which relates the elevation plane model coordinates to the reference 
coordinate system. The transformation of coordinate systems is 
examined in Appendix I. . . . 

Notice that the coordinate transformations described in 
Appendix I are very important for the analysis. This is especially 
true when the antenna under consideration is mounted off the 
center line of the fuselage. Using these transformations not only 
can the electromagnetic field in the desired radiation direction in 
the reference coordinate system be determined through the two 
principal plane cylinder solutions, but the polarization problem 
associated with the two different solutions can also be easily 
resolved without any difficulty. In addition, these transformations 
are the foundation on which the volumetric pattern analysis is based 
A detailed discussion on an off-center line mounted antenna is 
given in Appendix II. 

To determine the radiation angle (e^e* <)>re) in the elevation 
plane model coordinate system which corresponds to the desired radia 
tion direction (0ro> 4>ro) in the reference coordinate system, 

Eqs. (80) and (83) are employed. The subscripts e and o refer to 
the elevation plane model and reference coordinate systems, 
respectively. By definition, the unit radial vector for a given 
direction ( 0 , ({>) is given by 


(64) R = sin© coscf) x + sin 0 sin(j> y + cose z 


Hence, the spherical coordinates (0y.p> i)>^g) of the radiation 
vector are determined by 


(65) 




-1 ^re 
^re 



tan 



+ ) 
re 


1/2 


Z 

re 


with the components and found through 
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( 66 ) 
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For the case where (|) = 0° and e = 90° as seen in Fig. 36, 
Eq. (66) reduces to 


(67) 


-Ve" 


■ 1 0 o' 


- sin 6^^ cos 

''re 


0 0 1 


sin sin 



o 

1 

o 


cos J 


Finally, the radiated field in the desired radiation direc- 
tion can be determined from the field found in the elevation plane 
model analysis with an appropriate polarization modification. Thus, 
the solution obtained can be incorporated into the roll plane 
cylinder solution to form a complete solution for the volumetric 
pattern. 

Let us now consider the polarization effect due to the 
coordinate transformation, which is vital in achieving a com- 
plete solution. Recall that, in the far zone, the electric field 
can be written as 


(68) E = Eg e + J 
or 

(69) E = E' 6 ' + E' <^' 
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Fig. 36— niustration of the spherical 
coordinate system. 


where the prime system denotes the new coordinate system as seen 
in Fig. 36. Since the components E 0 and El are found from the 
elevation plane solutions, the corresponding E 9 and in the 
desired coordinate system are yet to be determined. If, one 
takes the dot product of Eq. (69) with the unit vector 6 , the 
9 component of the electric field results, i.e. 

(70) , Eq = E- (e- • e) + E' ii’ . 0 ) 

Similarly, the ((> component can be found as 

(71) E^ = E’ (9' • i) + E' (V . i) 


r 
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Thus, the polarization of the radiated field in terms of the 
reference coordinate system can be properly determined. 

Before the roll and elevation plane model solutions can be 
combined together, one must be concerned with the alignment 
of the surface normal at the antenna location in the two 
analyses. For an antenna mounted centrally on the fuselage, the 
surface normal points vertically upwards in the roll plane model 
analysis; whereas, in the elevation plane model it points in a 
direction normal to the profile surface at the antenna location. 

If the normals do not align properly, the two solutions would not 
match to give a complete and continuous solution. Consequently, 
the roll plane elliptic cylinder needs a tilt such that the normal 
direction is common to both solutions. It is noted that a radial 
monopole is mounted parallel to the surface normal at the antenna 
location on the aircraft fuselage. 


To accomplish this alignment, the X^oll sxis of the roll 

plane model coordinate system is chosen to line up parallel with 

the normal to the surface at the antenna location on the elevation 

composite elliptic cylinder as shown in Fig. 34. The tilt angle 

0 ^.,. is found to be 
tilt 


tan 


-1 


(72) 6 


tilt 


= < 


tan 


2 

^e ^s ref 
“e X, ref 


X, ,ef 


if Zg is negative 


if Z^ is positive 


where X 5 and Z^ pgf are the coordinates of the antenna source 
relative to the reference coordinate system. The quantities ag, 
bg, and bg are the necessary parameters for the composite elliptic 
cylinder model for the elevation plane model analysis. Notice that 
the origin of the roll plane cylinders is located a distance a^ 
directly under the antenna. Again, the quantity a^ represents an 
axis of the roll plane cylinder. 

It is also noted that, the parameter a^ is not equal to 
Xs ref which is the height of the antenna on the elevation 
composite cyl inder model . This is due to the fact that the 
two elliptic cylinders used to approximate the fuselage profile 
and cross-section do not match to give the true representation 
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of a three-dimensional aircraft fuselage. Hence, to obtain the 
necessary roll plane model solution, a coordinate transformation 
s i mi 1 ar to that of the -el evation“p1ane=model-=ana'lysts'=descri bed- 
above is needed. The required rotation angle (Bq, <|)o) of the roll 
plane model coordinate system is found to be 


(73) 


'tilt 





for a centrally mounted antenna. Since the Stilt ^ogle is 
rotated about the y'-axis, Eqs. (80) and (85) shall be used. 

One last note for the completion of our solution is the 
phase reference problem. For both elevation and roll plane model 
analyses, the origin of the individual coordinate system is 
chosen to be the phase reference point. Hence, the result 
obtained through these two principal plane model analyses needs to 
be modified so that a common phase center can be utilized. For 
convenience, the origin of the reference coordinate system is 
selected as the phase reference point for both models. In this 
way, the solutions blend smoothly together to form a complete 
result for the three-dimensional geometry under consideration. 

For each of the two solutions just described, only two 
rays are considered. This is due to the fact that there are 
only two dominant rays which propagate around an electrically 
large elliptic cylinder. Recall that four dominant rays 
contribute to the radiation pattern in the shadow region for an 
antenna mounted on a prolate spheriod. To account for these 
four rays necessary for the radiation pattern in the shadow 
region, both elevation and roll plane model results are needed 
in the pattern calculation. 

In addition, the three-dimensional effects of wings and 
stabilizers are, also, considered everywhere such that the rays 
reflected and diffracted from these scatterers can be included. 

The contribution of wings and horizontal stabilizers are 
handled in the same way as done previously. 


The effect of the cockpit/radome section is determined 
using flat plates attached to the fuselage just as the wings. 

Since computation time in analyzing a real cockpit radome 
section is So great, it is apparent that a simple model is needed. 
Originally, a finite bent plate approximation of the cockpit and 
nose section was considered because of an experimental investigation 
of radome. This resulted in a truncated fuselage model as 
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discussed in Reference [31]. However, the radiation pattern cal- 
culated using a flat plate model compared very favorably with the 
bent plate model. Hence, the flat plate model was adopted not 
only for its simplicity in analysis but, also, for its reduced 
computer running time. Notice that the plate is mounted at 
arbitrary angles relative to the cylinder to simulate the cockpit/ 
radome section as seen in Fig. 34 (a). This is an improvement 
over the previous roll plane model where the wings are restricted 
to be horizontal. 

Similar considerations were made in studying the vertical 
stabilizer which resulted in the selection of a bent plate 
model. Thus, the scattering effects of the cockpit and vertical 
stabilizer can be taken into account very simply. The theoretical 
model as shown in Fig. 34 (a) is an illustration of the analytical 
simulation of an actual aircraft. With all the necessary mechanisms 
being completed, the radiation and scattering by the fuselage, wings, 
cockpit, and stabilizers can now be calculated. Finally, the 
volumetric pattern of a fuselage mounted antenna is obtained 
by summing the direct field from the source and the reflected 
and diffracted field from the various scatters. 

Before it can be used to determine the volumetric pattern 
for an airborne antenna, this newly developed solution must be 
tested and verified. The elevation pattern of an axial slot 
mounted on a prolate spheroid, thus, was calculated and is 
presented in Fig. 37. The comparison between the measured and 
calculated results is very persuasive. This particular problem 
is an effective test in that it illustrates the accuracy for 
the limiting case of a small nearly spherical object as opposed 
to the large cylindrical aircraft model. 

To illustrate the validity of the new solution, the principal 
plane radiation patterns of antennas mounted on the KC135 aircraft 
that was previously computed are now computed with the improved 
solution. With the cockpit/ radome and vertical stabilizer being 
taken into account, the new computer model, consisting of an 80" 
by 3.75" right semi-ellipse and an 8" by 3,75" left semi-ellipse, 
is shown in Fig. 38. The cockpit/radome section is simulated 
by a flat plate and the vertical stabilizer is approximated by a 
bent plate model for better simulation. The elevation plane 
patterns for a short monopole mounted forward and over the wings 
are shown in Figs. 39 (a) and (b), respectively. The patterns 
for a circumferential KA-band waveguide are shown in Figs. 40 (a) 
and (b). The results for an axial KA-band waveguid are depicted 
in Figs. 41 (a) and (b). It is observed that the comparison 
between the calculated and measured results is very favorable. 
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slot mounted on a 4a x 2 a prolate slot mounted on a 4A x 2A prolate 
spheroid using the surface of rev- spheroid using newly developed ■ 
olution approach. (Geometry illus- volumetric solution. (Geometry ; 
trated in Fig. 30). illustrated in Fig. 30). j 
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Fig. 38a— Computer simulated model for the fuselage profile of a KC-135 
aircraft (side view). 
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Fig. 38b— Computer simulated model for the cross-section 
(at antenna location) of a KC-135 aircraft 
(front view). 
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Fig. 38c— Computer simulated model for a KC-135 aircraft (top view). 
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Fig. 39a — Elevation plane pattern for a x/4 monopole 
mounted forward of the wings on a KC-135 
aircraft (with radome and vertical stabilizer 
included). 
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„Fig. 39b--Elevation plane pattern for a a/ 4 monopole 
mounted above the wings on a KC-135 aircraft 
(with radome and vertical stabilizer included). 
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Fig. 40a— Elevation plane pattern for a circumferential 
KA-band waveguide mounted forward of the wings 
on a KC-135 aircraft (wi th racfome and vertical 
stabilizer included). 
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Fig. 40b— Elevation plane pattern for a circumferential 
KA-band waveguide mounted above the wings on a 
KC-135 aircraft (with radome and vertical 
stabilizer included). 
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Fig. 41a— Elevation plane pattern for an axial KA-band 
waveguide mounted fonvard of the wings on a 
KC-1 35 aircraft (with radome and vertical 
stabilizer included). 
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F1g. 41b — Elevation plane pattern for an axial KA-band 
waveguide mounted above the wings on a KC-135 
aircraft (with radome and vertical stabilizer 
included). 
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The discrepancy displayed in aft sectors in Figs. 14 and 15 no 
longer exist in Figs. 39. and 40^ For^axia=l-slot-cases,— ^ _ 
the back lobe is picked up as compared with the previous results 
in section A. Notice that the nulls of the monopole for the 
theory and measurement do not align together as one can see in 
Figs. 39 (a) and (b). This can be attributed to some extent 
to the misalignment of the monopole and surface normal in the 
experimental work. 

Finally, the roll plane radiation patterns for the new 
model are calculated. In each case, the pattern is only slightly 
improved over previously calculated results as shown in 
Figs. 27 and 29. This is attributed to the fact that the 
surface rays contribution from the elevation plane model is 
negligible in the plane of interest. In any event, these re- 
sults demonstrate the validity and applicability of the new 
solution in determining the radiation patterns for general 
fuselage mounted antennas. 

To further demonstrate the versatility of this solution, 
the radiation pattern for an infinitesimal monopole mounted on 
the bottom of the fuselage of a Boeing 737 with a dielectric 
radome is calculated. The model for this case is shown in 
Fig. 42 where the radome effect is simulated by a flat plate 
mounted vertically upward. This model is adopted because of the 
similar consideration as discussed previously. Even though 
double diffractions are not included in the calculation, the 
resultant pattern is still in very good agreement with the 
measured result as shown in Fig. 43. This, again, illustrates 
the capability of the solution in predicting the radiation pattern 
of fuselage mounted airborne antennas. 
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Fig. 42— Computer simulated model for a A/4 monopole mounted at station 222 
on the bottom of the fuselage of a Boeing 737 aircraft (with radome/ 
nose section being modelled by a flat plateK 
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Fig. 43— Elevation plane pattern of a x/4 monopole mounted 
at stations 222 on the bottom of a Boeing 737 air- 
craft. (<)> = 0° at the left; (fi = 180° at the 
right). ■ 


97 


TAIL 


CHAPTER IV 

VOLUMETRIC PATTERNS OF AIRBORNE ANTENNAS 


The theoretical approach to analyze the complete volumetric 
radiation pattern of fuselage mounted antennas was developed in 
the last chapter. This theoretical solution was based on the two 
•principal plane analyses presented earlier and it has been veri- 
fied by applying it to various cases' previously calculated using 
principal plane analyses. The results show great improvement 
over previous calculations based on comparisons with 
experimental measurements. Hence, the new volumetric solution 
is now applied to analyze the volumetric patterns of fuselage 
mounted airborne antennas. 

The aircraft model of most interest in our study is the 
Boeing 737 aircraft in. that extensive experimental work is avail- 
able. The volumetric patterns of this aircraft model is examined 
in this chapter. 

The models used to simulate a Boeing 737 aircraft are shown 
in Fig. 44, where the roll plane cross-section of the fuselage 
is approximated by a 65.86" by 43.3" elliptic cylinder for an 
antenna mounted at Station 220 on the top of the fuselage. The 
elevation profile is modeled by a composite elliptic cylinder 
with a = 48.72", b = 308.56", and b' = 1307.04". Notice that the 
cockpit nose section and vertical stabilizer are approximated by 
finite flat and bent plates mounted obliquely on the fuselage 
which results in a simple model for our aircraft as shown in Fig. 
44. 

In Fig. 44(c), the top view of our model is shown to 
illustrate the finite three-dimensional effect of an aircraft. 

The dotted line indicates the width of the cross-sectional 
(roll plane) cylinder used in the calculations. The models 
of the wings, cockpit nose section, and stabilizers are also pre- 
sented. The geometry is taken directly from the three principal 
views of the aircraft scale model. As presented earlier, the 
coordinate systems used for both the elevation and roll plane 
geometries are depicted in Fig. 44. 

The radiation patterns for a a/ 4 monopole mounted at 
Station 220 above the cockpit on a Boeing 737 aircraft are, 
then, calculated using the model just described with the 
analysis presented in the last chapter. The three principal 
plane results are shown in Figs. 45 to 47 and found to be in 
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Fig. 44a— Computer simulated model for the fuselage profile 
of a Boeing 737 aircraft (side view). The antenna 
is located at station 220 on top of the fuselage. 
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Fig. 44b— Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at station 
220 on top of the fuselage. 
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Fig. 44c— Computer simulated model for a Boeing 737 aircraft 
(top view). The antenna is located at station 220 
on top of the fuselage. 
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Fig. 45— Elevation plane pattern of a A/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

((|) = 0° at the left; <P = 180° at the right). 
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Fig. 46--Roll plane pattern of a A/4 monopole mounted at 
station 220 on top of a Boeing 737 aircraft. 

(<j) = 90° at the left; 4> = 270° at the right). 
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Fig. 47— Azimuth plane pattern of a A/4 monopole mounted 
at station on top of a Boeing 737 aircraft. 

(0 = 92°). 
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very good agreement with measurements. The experimental work 
was performed with a great deal-of—patienGe-and-preG-i s-ion by-— 
the teGhnical staff at NASA (Hampton, Va.), using a 1/llth 
scale model of Boeing 737 aircraft. 

The coordinate system used for the experimental measurements 
is shown in Fig. 48, in which the z-axis is vertically upward. 

In order to calculate a radiation pattern in terms of this 
experimental coordinate system, a transformation of coordinates 
is necessary so that the corresponding radiation direction in 
the analytical aircraft reference coordinate system can be 
determined. This is accomplished by considering the z-axis of 
the experimental coordinate system as a radial vector in the 
reference coordinate system. By inputting its spherical coordin- 
ates (9Zref ■ ‘I’Zy-ef = 0°) 3S the rotation angles 0^ and (j)^ 

required in the computer program, the corresponding radiation 
direction, in terms of the reference coordinate system, is 
obtained* 

Z 



Fig. 48-- Illustration of the coordinate system used for 
experimental measurements. 
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The complete volumetric pattern for a x/4 monopole located 
at Station 220 on top of a Boeing 737 aircraft is calculated. 

The off -principal plane elevation patterns are shown in Figs. 

49 to 52. The azimuth plane patterns are shown in Figs. 53 to 
60. In each case, the calculated result compare very favorably 
with the measurement. It is noted that the measurement results 
have some asymmetry in the patterns. This could be attributed to 
misalignment of the monopole with respect to the surface normal 
or the movement of the model due to shifting weight during the 
measurement. 

Finally, a volumetric gain pattern is presented in Fig. 61. 
The various directive gain regions are indicated by a color code. 
The directive gain, by definition [54], is 


Q/ \ _ radiation intensity (e, ( ))) 

' ’ average radiated power 

_ 47t [radiation intensity (e, <j>)1 
total power radiated 


The total power radiated by the antenna is obtained by numerically 
integrating the far field radiation pattern. Using Eq. (74), the 
gain function is determined for the entire space. The three- 
dimensional gain pattern is, then, plotted in color as shown in 
Fig. 61. For example, the pink color indicates the region in 
space where the gain level is greater than or equal to 0 dB. 

In other words, in this region the radiation intensity of the 
antenna is greater than or equal to that of an isotropic point 
source. The yellow color indicates the region where the radia- 
tion intensity is greater than -3 dB but less than 0 dB in gain 
level. Consequently, the yellow and pink regions represent the 
region in space where the gain level of the antenna is equal to 
or above -3 dB in comparison with an isotropic point source. 

Similarly, the green, orange, and purple stand for -6 dB, -10 dB, 
and -15 dB levels in gain, respectively. The experimental results, 
taken at NASA (Hampton, Va.) are shown in Fig. 62. The compari- 
son between the theoretical and experimental results is seen to 
be very encouraging. 

The agreement shown in the gain patterns demonstrates, 
indeed, the capability of the analytical solution in predicting 
the radiation patterns of airborne antennas mounted on complex 
aircraft structures. Through a modern high-speed computer, this 
analytic solution not only predicts the radiation patterns of 
on-aircraft antennas accurately but also very efficiently. This 
theoretical calculation saves a tremendous amount of expense and 
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manpower that is involved in building a scale model , mounting 
antennas, and taking. the-measurementST~ This is-one-of the ““ 
definite advantages of analytical solutions offer over scale 
model measurements in the antenna design approach. In addition 
the theoretical analysis enables an antenna designer to choose 
and evaluate the antenna types and locations for the desired 
coverage for specific applications in the aircraft design stage 
The relocation of an antenna at a future date could also be 
made easier through this analytical process. 



NOSE 



F1g. 49~Elevation plane pattern of a x/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

( 4 ) = 10° at the left; <P = 190° at the right). 
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Fig. 50--Elevation plane pattern of a x/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

((}) = 20° at the left; <P = 200° at the right). 
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Pig. 51— Elevation plane pattern of a a/ 4 inonopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

((j) = 30° at the left; <t> = 210° at the right). 
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Fig. 52— Elevation plane pattern of a A/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 

(if> = 40° at the left; <l> = 220° at the right). 
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Fig. 53--Azimuth plane pattern of a monopole mounted at 
station 220 on top of a Boeing 737 aircraft. 

(0 = 50°) 
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Fig. 54--Azimuth plane pattern of a V4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft, 
(e = 60°) 
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Fig. 56— Azimuth plane pattern of a A/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft, 
(e = 80°) 


115 



Fig. 57— Azimuth plane pattern of a A/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft, 
(e = 90°) 
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Fig. 58— Azimuth plane pattern of a A/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 
(6 = 100 °) 
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Fig. 59— Azimuth plane pattern of a A/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft, 
(e = 110°) 
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Fig. 60— Azimuth plane pattern of a A/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 
(0 = 120 °) 
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Fig. 61— Calculated volumetric directive gain pattern of a X/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. 
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Fig. 62— Measured volumetric directive gain pattern of a X/4 monopole mounted 
at station 220 on top of a Boeing 737 aircraft. This result was • 
performed at NASA (Hampton, Va.). 


CHAPTER V 

MICROWAVE LANDING SYSTEM AIRBORNE ANTENNAS 


Despite today's rapid advancement in aviation technology, 
the approach and landing problem remains one of aviation's lead- 
ing causes of accidents and fatalities. In order to improve the 
landing guidance and to meet the demands of increased air traffic 
flow around airports, the Federal Aviation Administration, in 
cooperation with international participation, is now developing 
a microwave approach and landing system known as the Microwave 
Landing System (MLS). This system will replace the existing 
lower frequency Instrument Landing System (ILS) now in operation 
at many airports around the world. The MLS is intended to 
significantly improve the guidance accuracy for the approach and 
landing of modern high-speed aircraft, especially during times 
of poor visibility under severe weather conditions. 

One of the main obstacles in the development of the MLS 
has been the determination of the appropiate polarization 
(horizontal or vertical) of the radiated electromagnetic field 
for efficient system operation. Most of the previous polariza- 
tion discussions on this matter have been based on ground effects; 
however, the airborne antenna is an integral part of the system and 
may have a definite influence on the polarization discussions. 
Another problem in the design of a safe, reliable and accurate 
approach and landing system is the location of the antennas on 
the aircraft structure in order to achieve the desired radiation 
coverage. 

In the past, scale model measurements were used to design 
and evaluate the antennas for a specific application. This 
approach of airborne antenna design process requires a great 
deal of engineering time and money. With the theoretical solu- 
tions developed in this study to predict the radiation pattern of 
fuselage mounted antennas efficiently via a high-speed computer, 
it is natural to apply them to analyze the MLS airborne problem 
in terms of the antenna complexity necessary to achieve the 
desired coverage. This is both a theoretical and experimental study 
in that various theoretical results are compared with scale model 
measurements. The experimental study is taking place at NASA 
(Hampton, Va.), where extensive scale model measurements for 
various aircraft are being performed. 
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As in any study of this type, it is initiated using a 
simple model of a basic aircraft neglecting the radome. Once 
this result is verified by scale model measurements, the radome 
is added to both the analysis and scale models. Again, these 
results are compared to verify the theoretical solutions. 

Initially, the study is based on the antenna performance in 
terms of a desired principal plane patterns. Finally, the 
newly developed solution of Chapter IV is applied to determine 
the complete volumetric pattern for fuselage mounted antennas. 

This is intended to examine the three-dimensional coverage of 
the antennas of interest based on the required radiation diagram 
specified by an FAA work force. 

With the numerical solutions verified by measurements, 
they can, then, be used to design the airborne antennas necessary 
to give the desired coverage for the MLS application. This is 
because the theoretical solutions can provide antenna designs 
much more efficiently than the scale model measurement approach. 
Furthermore, a parametric study of various typical aircraft 
shapes can be made using the verified numerical aircraft antenna 
solutions. The results should provide the polarization limita- 
tions based on the airborne antenna complexity. Note that this 
is a study of fuselage mounted antennas and not antennas mounted 
behind the radome. 

Initially, an antenna location on the bottom of an aircraft 
fuselage behind the radome was proposed for the MLS application. 

Since the MLS is for the approach and landing purpose, this 
location appeared to be very attractive for the forward cover- 
age requirement. An application of the volumetric solutions to 
the Boeing 737 aircraft revealed good agreement between theoretical 
results and NASA scale model measurements as shown in Fig. 63. 

The 737 aircraft was modeled by a composite elliptic cylinder with 
a 64.46" by 1200" right elliptic cylinder and a 64.46" by 232.52" 
left elliptic cylinder. The radome on the 737 aircraft as shown 
in Fig. 63 (b) is located approximately 186.97" from the origin 
and the antenna (a short monopole) is located 3 wavelengths away 
from the radome at a frequency of 3.18 GHz. The experimental results 
shown in these figures were obtained using an 1/1 ith scale model and 
performed at NASA (Hampton, Va.). In this case, vertical polari- 
zation seemed to have better forward coverage than horizontal 
polarization. 

Besides the radiation pattern results, the phase, or more 
specifically the phase center of the antenna, was of great 
concern to the FAA in the MLS application. If the phase center 
drifts too far from the actual antenna location, it would cause 
the system to appear to locate the aircraft at the wrong position 
in space. Generating phase data using experimental facilities 
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Fig. 63a— The scale model (1/11) of Boeing 737 aircraft 
in Anechoic Chamber at NASA (Hampton, Va.). 
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Fig. 63b— Computer simulated model for a a/4 monopole mounted at station 222 
on the bottom of the fuselage of a Boeing 737 aircraft. 
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Fig. 63c— Elevation pattern of monopole on a 1/llth scale 
model of Boeing 737 aircraft. 
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Fig. 63d— Phase error plot of monopole on a 1/llth scale 
model of Boeing 737 aircraft. 
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Fig. 63e — Elevation pattern of circumferential slot on a 
1/11 th scale model of Boeing 737 aircraft. 
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Fig. 63f — Phase error plot of circumferential slot on 
1/llth scale model of Boeing 737 aircraft. 
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Fig. 63g--Elevation pattern of axial slot on a 1/11 th 
scale model of Boeing 737 aircraft. 
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Fig. 63h— Phase error plot of axial slot on a 1/llth scale model of 
Boeing 737 aircraft. 


is extremely difficult, time-consuming and expensive. Alterna- 
tively, once the numerical solutions are verified based on 
amplitude measurements as previously done, they can be applied 
to generate the desired phase data with little or no additional 
effort. 

The phase error results from the interaction between the 
fuselage structure and antennas. This is an indication of the 
degree of distortion or effect the fuselage structure has on 
the performance of the antenna. It is the phase of the radiated 
field when the phase center of the antenna system is taken as the 
phase reference center. It is defined by 

error phase = actual phase - 1 ine-of-sight phase 

line-of-sight phase = free space phase of antenna 

without aircraft present. 

The radiation pattern along with the error phase plot may be 
used to determine the polarization selection of the MLS system. 
Additional results for various aircraft such as Boeing 707, B-1 , 
and F15A are presented in Ref. [31]. 

As a result of the study based on the principal plane per- 
formance, it is observed that horizontal polarization does not 
give the broad coverage obtainable with vertical polarization 
if just a single fuselage mounted antenna is used in the airborne 
system. The problem is that horizontally polarized fields are 
shorted out by the fuselage which confines the energy to sectors 
above the antenna rather than the broad azimuthal coverage as 
desired. On the other hand, based on the desired pattern per- 
formance as specified by an FAA work force, it also is doubtful 
whether a single vertically polarized antenna will function 
properly. In order to examine completely the antenna's perform- 
ance, one must investigate the complete volumetric pattern. 

As far as the phase data are concerned, the acceptable phase 
error is as yet unknown. It will not be discussed in this study. 

Note that the previous results did not include the landing 
gear effect. However, since the MLS system is used during 
landing, the presence of the landing gear must be taken into 
consideration in the pattern prediction, both analytically and 
experimentally. Thus, the calculated and measured volumetric 
patterns for a x/4 monopole mounted at Station 222 on the bottom 
of a Boeing 737 aircraft with the landing gear door included 
are presented in Figs. 64 and 65, respectively. The severe 
distortion in the radiation pattern displayed in Figs. 64 and 
65 shows the critical effect of the landing gear on the antenna 
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Fig. 64 — Calculated volumetric directive gain pattern of a x/4 monopole 
mounted at station 222 on the bottom of a Boeing 737 aircraft 
(with landing gear). 
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(with landing gear). This result was performed at NASA 
(Hampton, Va.). 



performance. It is this pattern distortion which forces one to 
seek an alternate location for MLS antennas.. . 

The theoretical simulation of the landing gear door is made by 
mounting two flat plates on the fuselage, using the computer model 
depicted in Fig. 66 where a 64.46" by 233.52" and 64.46" by 1200.0" 
composite elliptic cylinder model is employed for the aircraft 
fuselage. Although only the landing gear door is included in our 
calculation, the theoretical results compare favorably with 
the scale model measurements. 

Presently, the MLS antenna is being located in the area 
right above the cockpit section since it provides the desired 
forward coverage without being distorted by the landing gear, 
wings, etc. For the Boeing 737 aircraft. Stations 220, 250, and 
305 on the top of the fuselage are proposed as test locations for 
the antenna installation. In order to determine the best loca- 
tion for the MLS application, the elevation plane radiation 
patterns for a x/4 monopole mounted at these three locations are 
first obtained. The elevation plane patterns are used in that 
it is the most critical MLS sector. Both theoretical and 
experimental results from this study are presented in Figs. 67 
to 69. The computer simulated model of the Boeing 737 aircraft 
is shown in Fig. 70 in which a 78" by 211.6" and 78" by 850" 
composite elliptic cylinder is used to approximate its fuselage 
profile. However, different roll plane elliptic cylinders are used 
to approximate the various cross-sections at the antenna locations. 
The results reveal good agreement between the theoretical predictions 
and scale model measurements. In addition, the results indicate 
that Station 220 appears to be the best choice since it gives 
the best forward coverage which satisfies the MLS performance 
requirement. As seen in Fig. 74 the 0 dB sector in the antenna 
coverage requirement is the most important for the landing 
application. However, due to structural problems that exist on 
the actual aircraft fuselage, the antenna of interest cannot 
be actually mounted at Station 220 but must be moved to Station 
250. In addition, it must be mounted 4" off the fuselage center- 
line. 

To further evaluate the space coverage performance of the 
MLS antenna at location Station 250, the complete volumetric 
pattern needs to be examined. Scale model measurements for 
both polarizations, then, were performed at NASA, (Hampton, Va). 

The three-dimensional gain patterns as described in the last 
chapter were, also, obtained. Since vertical polarization 
appeared to be more promising for MLS applications, a simple 
A/4 monopole was the only antenna studied in the remainder 
of this research. 
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Fig. 66a- -Computer simulated model for the fuselage profile of a Boeing 737 
aircraft (side view). The antenna is located at station 222 on 
the bottom of the fuselage (with landing gear and radome present). 
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Fig. 66b— Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at station 
222 on the bottom of the fuselage with landing 
gear and radome present). 
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Fig. 66c— Computer simulated model for a Boeing 737 aircraft (top view). 
The antenna is located at station 222 on the bottom of the 
fuselage (with landing gear and radome present). 
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Fig. 67— Elevation plane pattern of a A/4 monopole 
mounted at station 220 on top of a Boeing 
737 aircraft (<p = 0° at the left; i> = 180° 
at the right). 
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Fig. 68--Elevat1on plane pattern of a A/4 monopole mounted 
at station 250 (off center) on top of a Boeing 737 
aircraft. (<t> = 0° at the left; 'P = 180° at the 
right). 
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Fig. 69--Elevation plane pattern of a x/4 monopole mounted 
at station 305 on top of a Boeing 737 aircraft 
((p=0° at the left; (j)=180'’ at the'right) . 
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Fig. 70a— Computer simulated model for the fuselage profile of a Boeing 737 
aircraft (side view). The antennas are located at station 220, 
250 (off-center), and 305 on top of the fuselage, respectively. 
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Fig. 70b--Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at station 
220 on top of the fuselage. 
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Fig. 70c--Computer simulated model for a Boeing 737 aircraft (top view). The 
antenna is located at station 220 on top of the fuselage. 
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Fig. 70d — Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at station 
250 on top of the fuselage (off center). 
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Fig. 70e--Computer simulated model for a Boeing 737 aircraft (top view). The 

antenna is located at station 250 on top of the fuselage (off center). 
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Fig. 70f--Computer simulated model for the cross section 
{at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at station 
305 on top of the fuselage. 
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Fig. 70g— Computer simulated model for a Boeing 737 aircraft (top view). The 
antenna is located at station 305 on top of the fuselage. 


The theoretical patterns for this configuration was obtained 
using our new volumetric-solution.- Figure 71 illustrates- the - . 

computer aircraft model used in the numerical calculations. As 
shown in Fig. 71(c), a bent plate with an inside angle of 11.4° 
is used to simulate the vertical stabilizer; whereas, a flat plate 
model is used to approximate the cockpit/radome section. For the 
sake of efficiency, no horizontal stabilizers are included in this 
model. Figure 72 illustrates the calculated three-dimensional 
(volumetric) gain pattern and Fig. 73 presents the results of 
the scale model measurements. The comparison between the calcu- 
lated and measured results is very encouraging. This indicates 
that the volumetric pattern solution developed in Chapter III 
not only can predict the complicated three-dimensional radiation 
pattern for antennas mounted on the center-line of the fuselage 
but, also, for antennas mounted off the center-line. 

Based on the MLS antenna coverage requirement as shown in 
Fig. 74, the radiation pattern of a single vertical polarization 
monopole, as seen in Fig. 73, is not quite adequate. Thus, a 
tail mounted monopole is proposed to be mounted at Station 950 
on the bottom of the Boeing 737 aircraft fuselage for missed 
approach purposes. This location is chosen because the slope 
of the fuselage surface at Station 950 is similar to that at 
Station 250 on the top of the fuselage s;uch that the space 
coverage of the two antennas tends to have a better match. To 
examine its space performance, the volumetric pattern of a x/4 
monopole mounted at Station 950 on the bottom of a Boeing 737 
fuselage was obtained. The three principal plane radiation 
patterns are shown in Figs. 75, to 77. The computer model, in 
this case, is composed of an elevation plane cylinder (a 110.312" 

X 548.325" and 110.312" x 552.72" composite elliptic cylinder) and 
a roll plane cylinder (66.55" x 66.55" elliptic cylinder) with 
four flat plates attached to the cylinders as shown in Fig. 78. 

Two of the four flat plates are used to simulate the wings and the 
other two simulate the horizontal stabilizers. The engines are 
neglected in the present model for simplicity. Finally, the 
volumetric gain pattern is presented in Figs. 79 and 80. It is 
noted that the calculated results, even without the engine effect 
included, compares very favorably with the scale model measurements. 

It is noted that for antennas mounted on top of an aircraft 
fuselage, the effect of the engine is negligible so that it can 
be excluded in the theoretical calculations as shown earlier in 
Reference [1]. However, the engine may have some effect on the 
radiation patterns when an antenna is mounted on the bottom of 
the fuselage, especially when the antenna is close to the wings. 

The effect of engines can be included in this solution using 
finite elliptic cylinders as presented in Reference [55]. In 
fact, this inclusion is an anticipated future effort. 
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Fig. 71a— Computer simulated model for the fuselage profile of a Boeing 737 
aircraft (side view). The antenna is located at station 250 on 
top of the fuselage (off center line). 
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Fig. 71b— Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at station 
250 on top of the fuselage (off center line). 
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Fig. 71c--Computer simulated model for a Boeing 737 aircraft (top view). The 
antenna is located at station 250 on top of the fuselage (off center 
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at station 250 on top of a Boeing 737 aircraft (off center) with 
landing gear. This result was performed at NASA (Hampton, Va.). 
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Fig. 74— MLS antenna space coverage diagram. (Each level indicates the 
minimum required dB level as compared to an isotropic source). 
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Fig. 75— Elevation plane pattern of a A/4 monopole mounted 
at station 950 on the bottom of a Boeing 737 air- 
craft (<l> = 0° at the left; <t> = 180° at the right). 
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Fig. 76 — Roll plane pattern of a A/4 monopoTe mounted at 

station 950 on the bottom of a Boeing 737 aircraft. 
(<j> = 90° at the left; <P = 270° at the right). 
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Fig. 77 — Azimuth plane pattern of a x/4 monopole mounted at 
station 950 on the bottom of a Boeing 737 aircraft 
(0 = 92°). 
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Fig. 78a— Computer simulated model for the fuselage profile 
of a Boeing 737 aircraft (side view). The antenna 
is located at station 950 on the bottom of the 
fuselage. 
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Fig. 78b--Computer simulated model for the cross section 
(at antenna location) of a Boeing 737 aircraft 
(front view). The antenna is located at station 
950 on the bottom of the fuselage. 
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Fig. 78c— Computer simulated model for a Boeing 737 aircraft (top view). 
The antenna is located at station 950 on the bottom of the 
fuselage. 
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Fig. 80— Measured volumetric directive gain pattern of a X/4 monopole mounted 
at station 950 on bottom of fuselage of a Boeing 737 aircraft without 
landing gear. This result was ' performed at NASA (Hampton, Va.). 





CHAPTER VI 

SUMMARY AND CONCLUSIONS 


High frequency solutions for the three-dimensional volu- 
metric radiation patterns of general type on-aircraft antennas 
mounted on the fuselage has been the object of this research. 

It is well known that the scale model measurement approach of 
airborne antenna design requires a great deal of engineering 
time and money. On the other hand, the theoretical analysis of 
on-aircraft antennas using a general surface of revolution model 
as done in Reference [2] results in a numerical solution which 
is very complicated, time-consuming, and uneconomical. Thus, 
an efficient analytical model for the prediction of the antenna 
system performance in the presence of complex aircraft structures 
is needed. This solution should be able to calculate the 
volumetric patterns of airborne antennas in an accurate and 
efficient manner such that it can be used to determine the 
antenna location and optimum antenna design for a given application. 
Such a solution has been developed in this study. 

The basic approach applied here is to break up the aircraft 
into its simplest structural forms. Due to the high frequency 
requirements of airborne antennas, these structures may be 
analyzed using ray optics techniques with numerical values 
obtained using GTD. Once the scattering from these various 
structures is found, it is then adapted to the aircraft model 
simply by adjusting the field incident on the structural scatterer. 
The only limitation of the GTD solution is that the source and 
various scattering centers be separated by at least a wavelength. 

In some cases, even this requirement can be relaxed. 

A theoretical solution was developed in Chapter IV to 
analyze complicated three-dimensional volumetric radiation pat- 
terns for fuselage mounted airborne antennas. This solution 
utilizes roll and elevation plane model analyses developed earlier. 
The procedure to combine these two solutions into a complete 
solution that can handle the volumetric pattern is based on a 
previous study of antennas mounted on prolate spheriods. The 
belt concept of blending these two solutions together is the 
key to the success of the complete solution. The use of flat 
or bent plates to approximate an aircraft cockpit nose section 
and vertical stabilizer is also new and quite useful, not only 
because of its simplicty in analysis but also due to its 
practicality in computation. 
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A numerical procedure that can be used to model a practical 
aircraft fuselage has also been,_presented. The-three-di mens tonal 
fuselage in the numerical model is simulated by two composite 
elliptic cylinders; one approximates the aircraft profile and 
the other its cross-section. These ellipses are numerically 
obtained through a best-fit ellipse routine. The description of 
the wings, cockpit, and stabilizers are measured directly from 
the three principal views of the aircraft. 

To demonstrate the validity and applicability of the 
theoretical solutions, the radiation patterns in the principal 
planes for antennas mounted on various locations on a KC135 
aircraft and a Boeing 737 aircraft have been calculated and 
compared with measurements taken at NASA (Hampton, Va.). The 
results show good agreement between the analytical and experi- 
mental Work and obvious improvement over previous analyses. To 
illustrate the versatility of the new solution, the radiation 
patterns for antennas mounted on a space shuttle with the radome 
taken into account have been calculated. Again, the result is 
very encouraging. 

The off-principal plane radiation patterns of a x/4 monopole 
mounted above the cockpit of a Boeing 737 aircraft have been calcu 
lated. The patterns obtained were shown to compare well with 
measurements taken at NASA (Hampton, Va.), which verify the 
assumptions made in the theoretical solutions. Color volumetric 
pattern plots were also presented in this study. The different 
colors in the pattern indicate various gain levels of the antenna 
pattern in comparison with an isotropic point source. These 
patterns provide the space coverage diagram necessary for the 
antenna designer. The success in predicting the volumetric 
patterns demonstrates the validity and capability of the 
theoretical solution to handle fuselage mounted airborne 
antennas. 

Finally, the volumetric pattern analysis of airborne antennas 
for the MLS application was presented in Chapter V. For the 
proposed antenna locations, the theoretical results compared very 
well with the scale model measurements. This, again, illustrates 
the applicability of the improved analytical solution in pre- 
dicting the complicated three-dimensional volumetric patterns 
for general type on-aircraft antennas. The study presented in 
Chapter V also showed that the vertical polarization seems to have 
broader space coverage than horizontal polarization. This is 
useful information for MLS designers. Note that the effect of 
landing gear is also presented to illustrate the necessity of 
critical and careful consideration in the selection of antenna 
location in a complicated modern aircraft. 
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The solutions developed in this research provide a useful, 
accurate, economical, and efficient means for determining the 
location and design of airborne antennas based on their pattern 
performance. For example, the program developed in this study has 
now been delivered to NASA (Hampton, Va.), and it typically 
runs a pattern in 30 seconds or less on a CDC 6600 digital 
computer. Consequently, this solution offers an excellent 
design tool which can be used to determine antenna location and 
designs for airborne applications either at the design stage of 
an aircraft or as a retrofit. 

Besides the radiation patterns, this numerical solution 
also provides phase data which may be of concern in the design 
of airborne antenna systems. Generating phase data using ex- 
perimental facilities is extremely difficult, time consuming 
and expensive. Since the desired phase data can be obtained 
with little or no additional effort, this is an additional 
advantage that the analytic solution offers over that of scale 
model measurements. Note that the accuracy of the computed 
phase data can be confirmed by the good agreement between com- 
puted and experimental amplitude data. 

It should be noted that for the antenna locations considered 
in the present study, the effect of the engines is negligible and 
they have been ignored in the model. However, the engines appear 
to have some effect on the radiation patterns when an antenna is 
mounted on the bottom of the fuselage. When necessary, the engine 
effect can be included in the solution by using a finite elliptic 
cylinder, as [>resented in Reference [55]. 
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APPENDIX I 

COORDINATE SYSTEMS TRANSFORMATION 


Consider that a vector R (x, y, z) in an xyz coordinate 
system as shown in Fig. 81 (a). A second coordinate system, 
(x'y'z') is also depicted, which is formed by first rotating 
the x-axis of the original coordinate system through an 
angle of <l)o about z-axis and then, the z-axis through an angle 
of 00 about the rotated x'-axis_^ Recall that in a vector space, 
the coordinates x-j of a vector X representing a point in 
n-dimensional spa^e can_also be regarded as the coefficients of 
the unit vectors e-j if X is represented as a sum of multiples 
of the unit vectors 


(75) X = 


- x^ e^ 


n n 


These unit vectors form a vector basis which span the n-dimensional 
vectoir space. In the three-dimensional geometry, where n = 3, 
let {x, y, z} be ^ basis of the vector space under consideration. 
Thus, the vector R can be written in terms of these unit vectors 
as 


(76) R=xx + yy+zz 


Let {x' , y'^ z'} be another basis of the same space such that 
the vector R is expressed as 


(77) R=x'x'+y'y'+z'z' 

Xhen, by a linear transformation, the coordinates of the vector 
R relative to the {x, y, z} b§sis^can^be determined from the 
coordinates relative to the {x‘, y', z'} basis through a matrix 
relation as 


(78) X = PX' 
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Fig. 81a— Illustration of coordinate system rotation. 

(Rotate <p angle about z axis and then rotate 
0Q angle about x' axis). 



Fig, 81b — Illustration of coordinate system rotation. 

(Rotate ({)q angle about z axis and then rotate 
6 angle about y' axis). 



Where X and X' are column vectors given by 




X 


■ x'“ 

(79) 

X = 

y 

and X ' = 

y‘ 



z 


z' _ 


and P^is a 3x3 transformation matrix from the {x, y, z) basis to 
{x* , y' , z'} basis. By multiplying^Eq. (78) by P-1 , the inverse 
of P, the coordinates of R" in^{x^, y' , z'} basis can be found in 
terms of the coordinates in {x, y, z} basis, i.e., 

(80) X' = P'^X ' 


Eqs. (78) and (80) describe the relation on which the coordinates 
of a vector in one coordinate system can be determined from its 
coordinates relative to another coordinate system. 


Since the unit vectors x' , y ' , and z' are related to the unit 
vectors, x, 9, and 5 in the following manner, i.e., 

^ ■ , A /S 

X* = COS <t>Q X + sin <()q y 

(81) y' = -cos 0Q sin <t>^x + cos cos (jj^y + sin e^z 

^ ^ A A 

z' = sin 9g sin<})Q x -sin 9^ cos y + cos 9^ z 

the transition matrix Px, which is defined as the transpose of 
the above matrix of coefficients, can be obtained as 


(82) P^ 


cos <|)„ -cos 9„ sin ((> sin 9„ sin <(>„■ 

0 0 0 0 0 

sin <()„ cos 9^ cos <(>„ -sin 9„ cos i 

0 0 0 0 0 

0 sin 9 - cos e„ 

o 0 J 


where subscript x indicates the 9 angle rotation is about the 
rotated x'-axis. By inverting the matrix Px, the inverse matrix 
Py"^ of P^ can be easily determined and is given by 

A X 
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0 



cos 

4> 

^0 

sin 

<^0 



-cos 

0 sin d> 

0 0 

cos 

% 

cos 

o 

-e- 

sin 

0^ sin <t> 

0 0 

-sin 

% 

cos 

^0 


sin 

0 

cos 0Q 


Similarly, if a new coordinate system, x'y'z' system, is set up 
by first rotating the x-axis of an angle of <f>o about z-axis and 
then the z-axis an angle of 6o about the y'-axis in a counter- 
clockwise sense as seen in Fig. 81 (b), the transformation matrices 
P and P "i can be determined as 

J J' 


(84) 



cos 

0 

0 

cos 

<t> 

-sin 

sin 

0 

0 

cos 

t 

O 

-e- 


cos 

% 

sin 

4 > 

^0 

cos (p^ 

sin 

®0 

sin 

^0 


-sin 

% 



0 

cos 

CD 

O 




and 


(85) 



cos 

% 

cos 

^0 

cos 


sin 

(|> 

^0 

-sin 0^ 

= 

-sin 

^0 



cos 




0 


_ sin 

0 

0 

cos 

<p 

0 

sin 

®0 

sin 

<f> 

^0 

cos 0Q 


where the subscript y indicates the rotation of 0o angle is about 
the rotated y'-axis in distinction with that of the x'-axis in 
Fig. 81 (a). 
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APPENDIX II 

GEOMETRICAL CONFIGURATION FOR OFF-CENTER LINE 
MOUNTED ANTENNAS 


Consider a source which is located at some angle off the 
center line of an aircraft fuselage as shown in Fig. 82. This 
source can be an infinitesimal monopole or arbitrarily oriented 
slot element. As discussed in Chapter III, the surface geometry 
nearest to the antenna has the most significant effect on the 
radiation patterns. For this case, the longitudinal surface 
curvature in the plane which is parallel to both the Z^ef axis 
and the surface normal at the source location, and the trans- 
verse surface curvature in the other plane which is orthogonal 
to the longitudinal plane, as shown in Fig. 82, play a dominant 
role in predicting the radiation patterns of airborne antenna. 
However, the transverse surface profile which cuts obliquely 
across the fuselage is not easily obtained from a scale model 
drawing of an aircraft (usually, only the elevation profile and 
cross section are given). Thus, to obtain the necessary elliptic 
cylinders for the volumetric pattern analysis, the aircraft 
cross section which cuts orthogonally through the source 
location is used to approximate the required transverse surface 
profile cylinder. This approximation is reasonable since the 
transverse surface profile of an actual aircraft does not 
change drastically for a small angle deviation. For the other 
elliptic cylinder which approximates the longitudinal surface 
profile at the source location, the following approach can be 
employed. 

First, let us approximate the elevation profile of the aircraft 
of interest by a composite elliptic cylinder as done previously 
in Chapter III. Its necessary parameters are given by Bq, semi- 
major axis, bo and bo', semi-minor axis for left and right half 
elliptic cylinder respectively, as shown in Fig. 83. Recall that 
the source location in our volumetric analysis is defined by 
(<l>s>Zs), as seen in Figs. 83 and 84, relative to the aircraft 
fuselage cross section where the source is located. Since the 
position of the source is given in terms of an elliptic cylinder 
approximating the fuselage cross section at the source location, 
its Cartesian coordinates can be easily determined as 
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Fig. 82— Geometri cal configuration of the longitudinal and transverse 
surface profile at the antenna location. 



Fig. 83— Front view of the roll plane (transverse) cylinder 
showing the antenna location and its surface normal . 
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Fig. 84 — Geometrical configuration of the antenna location. 


(86) 


where 


"s' ' Os h 


^s' = Os *s 


Os “ 


af bf 




2 2.2 
cos (f)^ + sin (j>^ 


The quantities a^ and are the semi -major and semi -mi nor axis 
of the cross sectional elliptic cylinder. In terms of the 
elliptic cylinder coordinate system as discussed in Chapter V, 


Xs and y^' are given by 


(87) 


where 


*5 ' V '^0= ''s 


bj; sin V. 


"s ' b; ♦s 


From Fig. 84, the distance Xeo between points P and Q is 
found via the ellipse equation as 
r 

a. 


(88) X 


^ b2 

N 0 


Zs if Zg is negative 


eo 


0 rT2 X 

'J^o ■ ^s positive 


where ao, bo» and bg have been defined earlier. Since the origin of 
the cross sectional cylinder is defined at a^ from point P along 
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line PQ, the source position in terms of the reference coordiante 
system can be easily determined from Figs. 83 and 84 as 


X 


s 



a^) + a^ cos Vg 


(89) 


= bf sin Vj 


Once the source location is determined in the reference 
coordinate system, the actual composite elliptic cylinder to 
simulate the longitudinal surface profile at the antenna 
location is given as 


(90) 


= < 


2x2 

^ ^s 

bn - 
0 s 


2 2 
X. + yt 
s •'s 


.2 


b„ = b. 


if is positive 


if is negative 



As described earlier, the elliptic cylinder which simulates the 
fuselage cross section at source location is used to approximate 
the transverse surface profile. Consequently, the necessary axes 
of the transverse elliptic cylinder are given by 



Note that the transverse surface profile discussed previously 
lies in a plane which is parallel to the surface normal at the 
source location. Hence, the transverse or roll plane elliptic 
cylinder needs a tilt such that the cylinder surface normal 
aligns with the original surface normal as seen in Fig. 85. 
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model relative to the reference coordinate system for 
center line mounted antennas. 


The required tilt or rotation ( 0 ^ tilt’ tilt^ given by 


r 


(92) 


®r tilt ’ 


tan — — = if z is positive 

2 f~2 2 ^ 


tan 


1 z 
■1 e s 


, 2 | 2 ^ 2 

'e 


if Zg is negative 


“^r tilt ■ ° 


with respect to the reference coordinate system. It is noted that 
the rotation described here is about y^gf axis and not axis. 

Notice also that the coordinate system of the elevation plane 
(longitudinal) model is rotated an angle ((|)n) which is the angle 
between the surface normal and Xref axis as shown in Fig. 83 , 

relative to the reference coordinate system. In other words, the 
elevation model coordinate system shown in Fig. 86 is obtained 
through a rotation <|)tilt) o'*" reference coordinate system 

as discussed in Appenaix I. The rotation angles are given by 


(93) 


where 


®e tilt 


90° 


"^e tilt ‘*’n 



tan 


a^ sin v^ 
f s 

bo cos V 
f s 


with respect to the reference coordinate system. The geometry of 
the coordinates is shown in Fig. 86. 


With the necessary elliptic cylinders being determined, the 
procedure to combine these two models to simulate the actual 
aircraft is the same as discussed in Chapter III. Thus, the 
volumetric patterns for off-center fuselage mounted airborne 
antennas can be determined. 
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model relative to the reference coordinate system for off center 
line mounted antennas. 
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